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CHAPTER  I 
INTRODUCTION 

In  the  design  of  the  cockpit,  an  important  constraint  is  the  feasibility  of  a 
safe  ejection  from  the  aircraft.  Thus,  in  the  consideration  of  new  designs  or 
design  modifications,  their  effect  on  the  safety  of  ejection  must  be  evaluated. 

For  example,  the  expected  posture  of  the  pilot,  pro-ejection  alignment,  and  the 
angle  of  ejection  all  influence  the  forces  sustained  by  the  spine.  Therefore  the 
designer  requires  information  as  to  which  combination  of  these  parameters  will 
minimize  the  possibility  of  injury.  Moreover,  it  is  useful  to  be  able  to  study 
particularly  detrimental  combinations  of  parameters,  so  that  the  ejection  system 
and  cockpit  environment  can  be  designed  to  minimize  injuries  even  under  extreme 
circum C ances. 

Experimental  work  on  ejection  has  consisted  largely  of  laboratory  tests  con- 
ducted on  human  volunteers  and  animals.  These  tests  are  expensive,  require 
considerable  preparation,  and  provide  only  limited  information  as  to  the  possibility 
of  injury.  Furthermore,  such  tests  are  not  often  possible  within  the  time  framework 
of  the  design  process.  Thus  the  design  of  ejection  devices  has  always  depended 
heavily  on  the  use  of  analytical  models.  The  most  well  known  of  these  is  the 
Dynamic  Response  Index  Model  (DR'i),  which  is  essentially  a one  degree  of  freedom 
lumped  mass,  spring  model  of  the  head-spine-torso.  It  has  been  extensively 
correlated  with  injury  data  and  provides  a useful  criterion  for  evaluating  the 
stfety  of  a crewman  exposed  to  axial  acceleration. 

For  situations  where  specifications  other  than  the  axial  response  must  be  con- 
sidered, simple  models  of  the  DRI  type  become  less  attractive,  because  each  para- 
meter, such  as  curvature  of  the  spine,  forward  tilt,  or  asymmetric  head  mass,  would 
have  to  be  correlated  with  injury  statistics.  For  these  applications,  general 
purpose  simulation  models  with  the  capability  of  treating  a wide  variety  of  environ- 


meats  appear  more  promising. 

Belytschko,  et  al.,  (1976,  AMRL-TR-76-10)  have  developed  three  dimensional 
models  of  the  head-neck-spine  which  treat  the  individual  anatomical  elements,  such 
as  vertebrae,  intervertebral  discs,  ligaments  and  ribs.  Because  the  geometry  and 
physical  properties  are  modeled,  specific  features  of  the  pilot's  anatomy  and  his 
equipment  can  be  treated.  Moreover,  the  models'  behavior  under  forces  in  any 
direction  and  of  any  intensity  can  be  studied,  so  that  various  environments  can 
be  evaluated. 

As  a consequence  of  its  generality,  two  aspects  of  the  model  need  attention: 
itscompletenesa  leads  to  considerable  computational  requirements,  and  it  requires 
validation  by  comparison  with  experimental  results. 

In  order  to  reduce  the  computational  requirements,  the  model  has  been 
modularised  so  that  .rU.es  components  can  be  simplified.  Four  models  have  evolved 
and  are  described  in  Chapter  II.  The  first  three  are  from  AMRL-TR-76-10,  although 
aome  data  has  been  modified  as  the  result  of  validation  studies.  The  fourth  model 
was  developed  completely  in  this  study.  It  is  a simplified  version  of  the  more 
general  model  of  the  spine,  yet  it  duplicates  most  aspects  of  the  more  complex 
models  under  vertical  excitation,  and  thus  proves  ideal  for  design  use. 

The  raw  output  of  these  models  are  force,  moment,  displacement,  velocity  and 
acceleration  time  histories  at  the  various  anatomical  elements.  For  the  purpose  of 
facilitating  the  interpretation  of  this  output,  an  injury  criterion  has  been  devel- 
oped, and  is  described  in  Chapter  IV.  This  injury  criterion  is  based  on  the 
stresses  in  the  cortical  bone  of  the  vertebral  body.  By  using  a statistical  dis- 
tribution about  the  mean  breaking  strength,  it  has  been  possible  to  develop  • post- 
processor package  that  gives  the  likelihood  of  failure  at  the  levels  of  the  spine 
for  each  simulation. 

The  validation  of  this  model  is  being  approached  in  two  ways.  First,  the 
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impedance  of  the  model  to  vertical  excitation. at  tb,e  pelvis  hae  been  compared  to  ,, 
the  experimental  measurements  pn  human  subjects.  To, obtain. the  model  impedance, 
a Vast  Fourier  Transform  technique  has  been  developed  which  enables  the  impedance 
of  any  linearized  model  to  be  found  from  a single  simulation.  This  procedure  is 
described  in  Chapter  III  and  comparisons  of  the  model  impedance  with  the  experimen- 
tal results  are  Dresanted.  By  adjusting  the  damping  cf  the  models  in  AMRL-TR-76-1C, 
it  has  been  possible  to  duplicate  experimental  impedance  curves  very  closely. 
Furthermore,  it  has  been  shown  that  the  major  peak  in  the  impedance  curve  at  approx- 
imately 5 He  arises  from  a combination  of  three  resonances:  the  seat-buttocks,  the 
flexural  response  oi  the  spine,  and  the  visceral  resonance. 

The  second  approach  which  has  been  taken  to  validate  the  modeling  procedure 
consists  of  studying  similar  models  "if  non-human  primates.  Since  experimental 
results,  particular'/  in  high  acceleration  environments,  can  more  readily  be 
obtained  for  animal  subjects,  this  will  enable  the  basic  modeling  procedures  to  be 

evaluated.  In  addition,  the  development  of  similar  models  for  non-human  primates 

/ 

will  hopefully  shed  soma  light  on  the  problems  of  scaling.  These  models  are 
described  in  Chapter  V.  It  is  of  interest  to  note  that  even  these  preliminary 
models  indicate  some  of  the  sources  of  difficulties  in  scaling  between  primates: 
for  example,  the  lowest  frequency  of  the  chimpanzee  model  is  dominated  core  by  the 
head-neck  mode  than  that  of  the  human  body  because  of  the  larger  relative  mass  of 


the  chimpanzee  head. 

It;  the  course  of  this  work,  the  problem  of  ejection  has  continuously  been 
studied  with  the  currently  available  models.  The  purpose  of  those  studies  has  been 
to  obtain  an  understanding  of  what  factors  play  an  important  role  in  the  response 
of  the  model  and  an  assessment  of  the  validity  of  the  model.  In  addition, 
concern  has  arisen  about  the  influence  of  the  behavior  of  the  spine  during  pre- 
ejection alignment  by  the  inertial  reel.  Several  studies  hsve  been  made  of  this 


U 


process  with  the  model  and  are  described  in  Chapter  VI.  These  studies  include  the 
effects  of; head  position,  the  restraint  system,  and  harness  positions.  In  Chapter 
VII,  ejection  simulations  are  reported.  These  include  results  from  models  of 


different  levels  of  complexity,  and  it  is  shovn  that  the  simplified  version  of  the 
spine  model  is  <;uite  representative  of  the  response  in  an  ejection  environment. 


...  ’ '■  • • . ~ * \ 
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CHAPTER  II  

SPINE  MODELS 

In  this  Chapter  the  geometrical,  stiffness  and  inertial  data  of  the  spine 
models  are  sunmarized.  Four  models  will  be  presented: 

1.  the  isolated,  ligamentous  spine  and  head  (ILS); 

2.  a two  load-path  model  consisting  of  the  isolated  ligamentous  spine,  the 
viscera,  and  head  (ILSV); 

3.  a complex  spine  model,  which  includes  the  spine,  rib  cage,  viscera  and 
head  (CSM); 

4.  a simplified  spine  model  which  uses  a minimal  number  of  elements  and 
masses  to  represent  the  spine,  viscera,  and  head  (SSM). 

The  first  three  models  have  previously  been  described  by  Belytschko,  et  al., 
(1976) , but  the  viscera  in  these  models  have  been  modified  extensively  on  the 
basis  of  impedance  studies  which  are  described  later. 

The  fourth  model,  SSM,  is  a new  model  which  is  aimed  at  providing  a simple 
design  tool.  This  model  has  only  8 nodes,  co*"oar.iJ  to  27  nodes  for  ILSV  and  42 
nodes  for  CSM,  respectively.  Running  time  for  an  100  msec  simulation  on  an  IBM 
370/158  is  about  10  sec,  which  is  of  the  order  of  1/40  of  the  computer  time  for 
ILSV  and  about  1/100  of  the  computer  time  for  CSM. 

The  geometries  of  the  models  are  defined  relative  to  a global  coordinate 
system;  its  origin  is  4 cm  directly  above  the  pelvic  mass  center.  The  positive  x 
axis  is  oriented  sideways  (to  the  left),  the  y axis  is  positive  towards  the  back, 
and  the  z axis  is  positive  vertically  upward  (towards  the  head).  The  x-z,  y-s, 
and  x-y  planes  correspond  to  the  frontal  plane,  the  sagittal  plane,  and  the 
horizontal  plane  respectively.  Many  of  the  procedures  used  to  determine  the  data 
were  discussed  in  detail  by  Belytschko,  et  al.,  (1976);  We  will  here  confine 
ourselves  to  alterations  of  the  data. 
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II.  1 isolated  Ligamentous  Spine 


The  Isolated  ligamentous  spine  model  consists  of  the  head,  the  cervical  and 
thoracolumbar  spines,  and  the  pelvis.  Figures  2.1  and  2.2  depict  the  frontal 
plane  and  sagittal  plane  views  respectively.  Tie  crosses  in  Fig.  2.2  represent 
primary  nodes,  which  coincide  with  the  mass  centers  of  the  rigid  bodies.  They 
are  generally  located  outside  the  vertebrae  because  they  represent  the  mass 
centers  of  entire  horizontal  torso  segments.  The  vertebrae  are  considered  rigidly 
embedded  in  their  respective  torso  segments.  This  representation  was  first 
conceived  by  Ome  and  Liu  (1970). 

The  geometry  of  each  vertebra  is  determined  by  the  position  of  13  points,  or 
secondary  nodes,  representing  the  centers  of  the  inferior  and  superior  end  plates; 
the  spinous  process  tip;  the  left  and  right  transverse  process  tips;  the  left  and 
right/inferior  - superior  articular  facet  points;  and  the  left  and  r ight/ inferior  - 
superior  ligaments  flava  points.  Adjacent  vertebrae  are  interconnected  through 
the  secondary  nodes  by  deformable  elements  representing  the  connective  tissues; 
ligaments  and  the  intervertebral  disc.  Spring  elements  with  resistance  only  in 
tension  are  used  to  represent  the  ligaments.  Although  data  defining  the  articular 
facet  planes  has  been  developed,  in  the  thoracolumbar  spine  the  facet  joints  ore 
represented  by  spring  elements. 

In  the  cervical  spine,  each  articular  facet  plane  is  represented  by  three 
points  and  the  joints  are  modeled  with  hydrodynamic  elements.  The  intervertebral 
discs  are  modeled  as  short  beams.  These  elements  are  described  in  detail  by 
Belytschfco,  at  al,,  (1976),  Tables  2.1,  2.2  and  2.3  list  the  geometries,  stiffness 
and  inertial  data  for  the  isolated  ligamentous  spine.  The  inertial  data  Is  based 
o©  Liu  and  Wiekstroca  (1973),  the  stiffness  data  primarily  on  Schultz,  et  a l , 
(1973),  except  that  the  loabar  disc  axial  stiffnesses  are  increased  to  account  for 
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Table  2.1  Isolated  ligamentous  spine  (ILS)  geometric  data  ss  reported 
by  Belytschko,  et  si.,  (1976). 


Vertebral 

Level 

Coordinates  of 

Inferior  End  Plate  Center 

y (cm)  2 (cm) 

Vertebral  Body 
Height 

(cm) 

Intervertebral 
Disc  Height* 

(cm) 

L5 

1.800 

2.020 

2.392 

1.859 

L4 

1.100 

5.700 

2.636 

1.354 

L3 

1.000 

9.550 

2.751 

1.223 

L2 

1.331 

13.450 

2.792 

1.173 

LI 

2.142 

17.150 

2.726 

0.996 

T12 

3.003 

20.590 

2.567 

0.822 

Til 

3.882 

23.680 

2.433 

0.645 

T10 

4.594 

26.500 

2.298 

0.477 

T9 

4.849 

29.240 

2.146 

0.460 

T8 

4.638 

31.830 

2.073 

0.459 

T7 

4.580 

34.300 

2.019 

0.404 

T6 

4.250 

36.610 

1.990 

0.314 

T5 

3.990 

38.850 

1.957 

0.266 

T4 

3.690 

41.000 

1.902 

0.214 

T3 

3.350 

43.150 

1.850 

0.274 

T2 

2.920 

45.260 

1.790 

0.306 

T1 

2.410 

47.440 

1.648 

0.448 

C7 

1.909 

49.420 

1.612 

0.394 

C6 

1.760 

51.448 

1.516 

0.434 

C5 

1.460 

53.516 

1.515 

0.576 

C4 

1.290 

55.439 

1.513 

0.417 

C3 

1.484 

57.332 

1.511 

0.398 

C2 

1 . 636 

59.239 

1.500 

0.408 

* Indicates  disc  below  vertebral  level. 


Table  2.2  Isolated  ligamentous  spine  (ILS)  Stiffness  Oats  from  Balytschko, 
et  al.,  (1976). 


Disc 

Level 

Axial 

Stiffness 

g 

Dyne/cm  x 10 

Torsional 

Stiffness 

9 

Dyne -cm  x 10 

Bending  Stiffness 
Dyne-cm  x 10® 
Sagittal  Frontal 

Plane  P7.ane 

Shear 

Deformation 

Parameter  4> 

S-L5 

1.47 

0.90 

0.70 

1.57 

7 

L5-L4 

1.87 

1.10 

0.80 

1.81 

9 

L4-L3 

- 2.00 

1.20 

0.90 

2.19 

9 

L3-L2 

2.00 

1.20 

0.90 

2.16 

13 

L2-L1 

2,13 

1.20 

0.90 

2.20 

14 

L1-T12 

1.80 

1.00 

0.90 

2.27 

16 

T12-T11 

1,50 

0.80 

1.00 

2.34 

30 

T11-T10 

1.50 

0.70 

1.20 

2.44 

41 

T10-T9 

1.50 

0.70 

1.10 

1.93 

72 

T9-T8 

1.50 

0.60 

1.10 

1.76 

78 

T8-T7 

1.50 

0.60 

1.00 

1.59 

66 

T7-T6 

1.80 

0.60 

1.00 

1.62 

82 

T6-T5 

1,90 

0.60 

1.00 

1.68 

•177 

T5-T4 

2.10 

0,60 

1,00 

1.81 

158 

T4-T3 

1.50 

0.40 

0.60 

1,22 

128 

T3-T2 

1.20 

0.30 

0.40 

1.04 

48 

72-T1 

0.70 

0.20 

0.20 

0,63 

34 

Tl-C? 

1.10 

0.18 

0.20 

0,62 

45 

C7-C6  V 

2.84 

0.11 

0.11 

0.25 

45 

C6-C5 

1.01 

0.07 

0.06 

0.10 

' 45 

C5-C4 

0.60 

0.04 

0.04 

0.06 

45 

C4-C3 

. 0.67 

0.03 

0,04 

0.05 

45 

C3-C2 

0.70 

0.05 

0.04 

0-05 

45 

C2-C1 

0.60 

0.04 

0.04 

0.10 

45 

Kaek*-l) 

0.14 

0.52 

0,6  9 

0,69 

2 

7 

MI  ligament  and  facet  axial  stiffness  1.5  x 10  dy^e/cm. 


'^Single  neck  eUsarot  which  c*n  bs  substituted  for  the  cervical  spine. 
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the  body-weight  preload  as  described  by  Belytschko,  et  el.,  (1976). 

In  many  studies  reported  herein  the  cervical  region  of  the  isolated  ligamen- 
tous spine  has  been  replaced  by  a single  beam  element*  the  stiffness  of  which  is 
given  in  Table  2.2.  The  isolated  ligamentous  spine  with  the  cervical  region 
represented  by  a single  beam  element  is  the  spine  model  which  will  be  referred  to 
as  the  ILS. 

The  ILS  model  allows  for  only  one  path  of  force  transmission  along  the  torso, 
namely  the  spine.  Some  investigators  (e.g.  Weiss  and  Mohr  (1967))  have  reported 
the  observation  of  a longitudinal  wave  through  the  abdomen  during  experiments  in 
which  subjects  seated  in  chairs  were  dropped  onto  yielding  as  well  as  non-yielding 
bases.  Because  the  spine  is  several  orders  of  magnitude  stiffer  than  the  viscera, 
it  is  the  primary  path  of  force  transmission  from  the  seat  to  the  upper  torso  and 
head.  However,  the  viscera  transmits  some  of  the  force  to  the  mid  and  upper 
torso  through  the  diaphragm  and  ribs.  It  is  these  observations  which  led  to  the 
development  of  the  two  visceral  models  which  will  be  discussed  shortly. 

Another  shortcoming  of  ILS  is  the  low  flexural  stiffness  of  the  spine  model. 
Lucas  and  Brassier  (1961)  determined  frontal  plane  buckling  loads  of 
2 x 10^  to  10  x 106  dynes  for  an  isolated  ligamentous  spine  constrained  against 
sagittal  plane  motion.  The  sagittal  and  frontal  plane  bending  stiffnesses  are 
of  the  same  order,  so  the  sagittal  plane  curvature  would  tend  tc  lower  the  sagittal 
plane  buckling  loads.  Since  the  spine  model  stiffnesses  are  based  on  the  material 
constants  of  the  intervertebral  discs,  it  will  also  buckle  under  these  loads. 

Tho  ejection  environment  produces  internal  forces  several  orders  of  magnitude 
higher  chan  the  buckling  loads.  To  increase  the  stability  of  the  spine  when 
subjected  to  these  accelerations  it  is  necessary  to  increase  the  beading  stiffness 
of  the  model.  This  is  accomplished  by  incorporating  representations  of  the  rib 
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cage and  viscera  into  the  spine  model. 

The  effects  of  the  viscera  and  the  rib  cage  are  approximated  by  a second 
column  of  beam  elements  interconnecting  the  primary  nodes  from  the  pelvis  to  Tl. 

The  properties  of  these  beam  elements  were  chosen  so  that  they  offer  resiatance 
only  in  cases  of  large  relative  rotation  between  the  ends  of  an  element  by  using 
a cubic  moment -curvature  relation.  As  will  be  Seen  subsequently,  this  secondary 
column  improved  Che  flexural  response  of  the  spine,  but  deformation  in  the  lumbar 
region  remained  quite  severe. 

II. 2 Isolat3d  Ligamentous  Spine  with  Viscera 

The  model  of  the  isolated  ligamentous  spine  with  a visceral  representation  will 
be  referred  to  as  the  ILSV  model.  In  ILSV,  vertebral  levels  Til-Si  are  separated 
into  visceral  and  spinal  masses.  The  subdivision  of  the  torso  cross-sections  la 
based  on  Eycleshymer  and  Schoetnaker  (1970) . A typical  section  and  the  areas 
ascribed  to  the  viscera  and  spine  are  shown  in  Fig.  2.3.  The  spinal  and  visceral 
translational  and  rotational  masses  at  each  vertebral  level  were  computed  from 

3 

these  areas  using  a uniform  density  of  1 gm/cm  , and  are  given  to  Table  2.4.  The 
total  masses  compare  well  with  the  data  of  Liu  and  Uickstrom  (1973). 

It  was  assumed  that  the  only  mode  of  force  transmission  through  the  viscera 
is  axial.  Hence  it  was  only  necessary  to  determine  the  axial  elastic  and  viscous 
constants  for  each  element  interconnecting  the  visceral  masses.  As  approximate 
fundamental  frequency  was  computed  for  the  viscera  using  the  results  of  Coermann, 
ct  al.,  (I960),  who  observed  a peak  in  abdominal  wall  displacements  when  the 
driving  frequency  was  between  3 and  4 Hr,  Asserting  that  displacements  of  the 
viscera  are  coupled  to  the  abdominal  wall  displacements,  with  peaks  occurring  at 
the  seats  frequencies,  we  esn  determine  the  fundamental  frequency  of  the  viscera. 

For  a damping  ratio  less  than  0.707,  the  frequency  ratio  which  defines  the  maximum 
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point  of  the  displacement  versus  excitation  frequency  curve  is 


where 

@/2t?  s fundamental  natural  frequency  (Hr) 

(u/2n  * excitation  frequency  (Hz)  ; 
cp  ® visceral  damping  ratio. 

From  the  results  of  Weiss  and  Mohr  (1967),  an  approximate  value  of  0.5  was 
determined  for  the  visceral  damping  ratio.  Using  a value  of  3.5  Hz  for  uj/2tt, 

Eq.  (2.1)  results  in  $/2tt  approximately  equal  to  5 Hz. 

If  we  treat  the  viscera  - torse  wall  system  as  a one-dimensional  rod  and 
assume  that  the  viscera  are  essentially  fixed  at  the  pelvis  and  the  diaphragm, 
the  fundamental  natural  frequency  is  given  by 

P/2tt  = 2^  (2.2) 

where  L » length  of  visceral  "stack" 

c «*  longitudinal  wave  speed  through  viscera. 

The  longitudinal  wave  speed  is  given  by 

c aJ~  (2 .3) 

V p 

where  B = effective  bulk  modulus  of  viscera  - torso  wail  system 

3 

p ® visceral  density  ■ 1 gm/cm  . 

Qombiniug  Eqa.  (2.2)  and  (2.3)  we  have 

B •»  (L0/fl)2  . (2.4) 

HithL  equal  to  31.64  cm  (the  distance  frcsa  the  pelvis  to  T10)  and  0 * lOn, 

5 2 

Eq.  (2.4)  yields  a value  of  1 x 10  dynes/ca  for  die  effective  bulk  modulus  of 
the  viscera  -^abdominal  wall  system. 

Once  we  have  determined  B,  the  stiffness  of  a visceral  element  is  given  by 


24 


where  A is  the  average  cross-sectional  area  of  the  element. 

The  XLSV  model  is  shown  in  Fig.  2.4.  At  vertebral  level  L5  to  Til,  the  visceral 
and  vertebral  mass  are  interconnected  by  horizontal  spring  elements  representing  the 
interaction  between  the  spine  and  the  viscera,  end  vertical  connecting  elements  with 
linear  axial  stiffness  and  cubic  bending  stiffness  were  used  to  connect  adjacent 
visceral  levels.  Since  there  is  no  data  for  the  viscere-spine  interactiwc,  the 
stiffnesses  of  these  elements  were  token  as  equal  for  all  levels  and  identified  by 
matching  whole  body  impedance  as  described  in  Chapter  III. 


Table  2.4  gives  the  mass  data  for  the  pelvis  through  T10  of  the  ILSV. 

Table  2.5  lists  the  stiffnesses  of  the  visceral  elements.  The  masses  for  the 
remainder  of  the  model  and  the  stiffnesses  of  the  intervertebral  discs  are  given 
in  Tables  2.2  and  2.3. 

XI. 3 Ligamentous  Spine  *rith  Ribs  and  Hydrodynamic  Viscera 

The  complex  spine  model,  which  will  be  referred  to  as  the  CSM,  includes  a 
detailed  representation  of  the  rib  cage,  consisting  of  rib  pairs  1 - IQ  and  the 
sternum.  Each  rib  and  the  sternum  are  modeled  as  a rigid  body.  The  ribs  interact 
with  vertebrae  Tl-TlO  through  the  costovertebral  and  costotransverse  Joints,  which 
are  modeled  as  spring  and  beam  elements  respectively.  Rib  pairs  1-7  are  connected 
to  the  sternum  through  the  costo-sternal  joints  which  are  modeled  as  beam  elements. 

Rib  pairs  8-10  arc-  connected  to  the  adjacent  ribs  by  spring  elements  representing 
the  behavior  of  the  in ter chondral  cartilage.  The  actions  of  the  intercostal  tissues 
are  al«o  represented  by  spring  elements. 

The  geometries  and  deformation  characteristics  of  the  rib  cage,  as  reported  by 
Belytschko,  et  el.,  (1976),  were  based  on  the  measurements  of  Schultz,  et  ai.,  (1974), 


Table  2.S.  Isolated  ligamentous  spine  with  viscera 
(IlSV)i  viscera  stiffnesses. , 


ttfy  slpfe 


*«*  •*;  Viscera  l^-.  ? ».? « 

Element 

> / ' * *s«f?b*w % 

(cm2) 

(cm) 

'Stiffness  v'  J 

(dyne/cm  x 10  ) ; i^A^vy  , 

■ Pelvia-Sl  •-«?  : 

191.56 **■ 

V 3.59 

• \‘[4 

' 0.92 

-v  ~2 6 3 ;7  7-  - - V-' 

• ;3.60-  — 

1 txii  TJ  . ■ ,(■'«  A 4 >'k’  {55  A*-  • - . »-  • 

■ - 1.29  * *-— ■ ■ ~ - - — -■  - S 

L5-L4 

370.80  -• 

3.83 

, Lv168  '■  j 

L4-L3 

- 371,35 

3.93 

1.64  . j 

L3-L2 

425.65 

3.87 

1.91  | 

L2-L1 

468.90 

3.65 

2.24  ' ' " -1 

•31 

LI-T12 

472.10 

3.37 

2.42  J 

T12-T11 

459.20 

3.03 

2,62  ' | 

«... 

Tll-TlO 

459.20 

2.79 

2.86 

All  viscera-spine  interconnecting  element  stiffnesses  » 1.0  x 10^  dyne/cm 


and  the  work  of  Andriacchi,  et  al.y  <1974)  , respectively.  - 


The  abdominal  viscera  are  modeled  by  hydrodynamic  elements  stacked  in  aeries 

5 2 

between  the  pelvis  and  rib  pair  10.  An  effective  bulk  modulus  of  lx  10  dynes /cm 
was  used  for  the  viscera.  Figures  2.5  and  2.6  depietthef rental  and  sagittal  plane 
views  of  the  ligamentous  spine  with  ribs  and  hydrodynamio  visceia.  Tables  2.6  and 
‘2.7  present  the  inertial  and  stiffness  data  as  reported  by  Belytschko,  et  al.,  (1976). 
The  stiffnesses  of  the  spine  elements  are  identical  to  that  given  in  Tabl&_.2.2< 

II.' 4 Simplified  Spine  Model  . 

Although  the  ILSV  is  the  model  of  primary  interest  for  the  impedance  study  in 
Chapter  III,  it  is  felt  that  even  this  model  is  computationally  too  time  consuming 
to  serve  as  a design  tool.  Therefore  an  equivalent  simplified  model  of  the  ILSV 
was  constructed.  This  simplified  spine  model  (Fig,  2.7),  will  here  be  referred  to 
as  SSH.  In  this  model,  the  thoracolumbar  spine  is  modeled  with  three  beam  elements, 

Tl  - TlO,  T10  - L3,  and  L3  - SI.  The  longitudinal  actions  of  the  viscera  are 
modeled  with  four  sping  elements  and  3 masses  between  the  pelvis  and  TlO.  The 
interaction  between  the  viscera  and  the  spine  is  modeled  with  a spring  element  at 
each  visceral  mass.  As  with  the  ILSV,  the  viscera/spine  interaction  spring 
stiffnesses  were  adjusted  by  the  impedance  studies. 

Tables  2.8  and  2.9  list  the  SSH  mass  and  stiffness  data.  The  SSH  translational 
mass  data  va9  computed  by  summing  the  appropriate  ILSV  translational  masses.  The  SSH 
rotational  mass  data  w?s  computed  using  the  parallel  axes  theorem.  The  SSM  stiffness 
data  was  determined  from  series  combinations  of  the  corresponding  ILSV  elements. 
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Table  2.6,  Complex  Spine  Model  (CSM>  inertial  data,  Belytacbko,  et  al,>  (1976) 


Vertebral 

level 

Maas  in 

g } - 

Grams  x 10 

Xxx  r>  A 

Gram-cm  x 10 

I— 

* yy  2 

Gram-cm 

V"  :•  5 

x 10* 

I— 

zz 

Gram-cm  x 

IQ5 

\ • . 

j 

Sv'* 

Pelvis: 

16.200 

128.000 

20.000  ; 

19.300 

LS 

' 1.500 

2.783 

1.795  \ 

Vv  ■ 

2.382 

I- 

L4 

1.500 

2. 748 

1.704 

■ v . i .If 

2.291 

•*..  . • 

L3 

1.500  , 0 

2.809 

. -1.682  : 

Jit  ;'1  f *• ' 

2.280 

L2 

1.500 

2.840 " 

1.695 

v*i  * 

Ic/C..'. 

2.291 

. ’ .*>" 

H 

1.500 

2.740 

1.569 

“ 

2.212 

. 

T12 

1.556 

7.002 

1.309 

1.919 

■ ■$£  ■ 

Til 

1.453 

7.056 

1.230 

1.941 

T10 

1.202 

6.028 

1.129 

■ t 

1.648 

so 

T9 

1.267 

6.164 

1.230 

1.716 

T8 

1.176  - 

S.543 

1.208 

<•  *'■ 

1.670 

’"T- 

T7 

1.158 

5.347 

1.219 

•V- 

1.6S9 

T6 

1.043 

4.425 

1.162 

1.546 

T5 

1.025 

3. 383 

1.151 

, . . 

1.490 

x 

T4 

0.964 

3.138 

1.060 

1.354 

T3 

1.010  ~ • 

2.878 

1.174 

- ■■ 

1.422 

T2 

0.974 

2.00? 

1.029  : 

1.230 

T1 

1.209 

0.745 

0.S18 

1.716 

. ‘ “ f . 

C7-C2 

1.000 

0.700 

0.500 

1.500 

Head 

5.612 

44.786 

4,044 

..  ; ■ ■ 

3.385 

Ribs  T1 
to  T10 

0.074 

0.373  ; . 

0.074 

.V  ■■ 

0.074 

Lower 

Viscera 

1.500 

10.700  7 

0.S50 

1.000 

' 1* 

Upper 

1.S00 

10.700 

0.550 

1.000 

Upper 

Viscera 


Tflblc  2.7  Stiffness  data  for  deformable  elements  in  CSM  rife  cage,  Bely tschko,  et  al.,  (197^ 


Table  2.3.  Simplified  spine  model  (SSM)  inertial  data 

Global  Coordinates^  Translational  Rotational  Hasses 

(cm)  Mfl88  (grama -cm2  x 105) 


Level 

y 

z 

3 

(grams  x 10  ) I-- 

I-- 

yy  ... 

Spine  Inertial 

Data 

Pelvis 

0.0 

-4.000 

> 14.880 

13.08 

20.29 

19.38 

L3 

4.497 

10.940 

1.660 

0.40 

0.58 

0.28 

T10 

1.469 

27.640 

- i 6.846 

5.328 

8.560 

8.307 

T1 

0.857 

48.240 

5.591 

2.773 

6.438 

7.212 

Head 

1.660 

66.800 

5.612 

4.479 

4.044 

3.385 

Visceral 

Inertial  Data 

ZJ> 

-3.407 

3.190 

4.295 

1.272 

2.985 

3.158 

L3 

-5.321 

10,940 

X 368 

1.203 

3.484 

3.746 

Ll 

-3.515 

18.460 

5.134 

2.576 

5.598 

6.257 

^Global  coordinates  of  mass  centers.  All  global  x coordinates  * 0.0 
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Tfible  2.9.  Simplified  spine  model  (SSM)  stiffness  date 


Element 

Axial 

Stiffness 

dyne /cm  x 10^ 

Torsional 

Stiffness 

g 

dyne/cm  x 10 

Bending 

Stiffness 

g 

dyne/co  x 10 

Neck 

14.40 

5.17 

6.89 

T1-T10 

15.40 

0.49 

0.67 

T10-L3 

34.70 

1.87 

2.15 

L3-S1 

56.80 

3.50 

3.43 

T10-L1  V<X> 

0.88 

LI-L3  V 

1.03 

L3-L3  V 

0.83 

LS-Pelvis  V 

0.54 

^ V a viscera 

7 

All  viscera-spine  interconnecting  element  stiffnesses  = 1.00  x 10  dyne /cm 

CHAPTER  III 


IMPEDANCE 

III.l  Rationale  and  Overview 

Driving  point  impedance  can  be  broadly  defined  as  the  ratio  of  the  force  at 
the  driving  point  to  the  velocity  of  the  driving  point.  More  specifically, 
impedance  is  the  complex  ratio  of  the  transforms  of  the  force  and  velocity.  The 
quantity  of  interest  to  us  is  the  impedance  magnitude  or  the  modulus  of  impedance, 
-which  is  a function  of  the  excitation  frequency.  An  impedance^  plot  for  a linear 
system  gives  insight  int'  the  frequency  content  of  the  system,  and  the  effect  of 
the  damping,  stiffness  and  mass  associated  with  different  elements  of  the  system. 

If  a model  is  to  duplicate  the  dynamic  behavior  of  a physical  system,  it  must 
duplicate  the  impedance  of  that  system.  In  particular,  any  model  of  the  human  body 
must  have  the  same  impedance  characteristics  as  that  measured  for  the  human  body 
if  any  validity  is  to  be  ascribed  to  it. 

It  i3  interesting  that  in  spite  of  the  fact  that  extensive  impedance 
measurements  have  been  made  of  the  human  body  and  in  fact  constitute  the  largest 
available  set  of  experimental  data  for  dynamic  response,  the  impedance  of  models 
such  as  chose  of  Orne  aud  Liu  (1970), and  Prasad  and  King  (1976)  were  not  compared 
to  experimental  results. 

One  Reason  for  this  is  that  the  determination  o*  the  impedance  through  direct 
\ 

integration  is  a laborious  process,  since  the  model  has  to  be  excited  harmonically 
until  steady  state  response  is  obtained  at  each  frequency  that  the  impedaac  is 
desired.  In  this  investigation*  a Fast  Fourier  Transform  technique  va*  developed 

TvT  “ 

'“'We  will  use  impedance  to  mean  modulus  of  impedance  for  convenience. 
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through  which  the  impedance  over  a frequency  ra»g<?  of  interest  can  be  determined 
from  one  time  history,  so  that  a single  run  suffices  toesteblish  the  impedance  -> 
curve. 

' \ . X 

The  mathematical  background  of  this  method  and  some  basics  of  the  concept  of 
impedance  are  developed  in  the  next  Section.  We  then  review  the  experimental 
literature  on  impedance  of  human  subjects.  The  final  two  .Sections  give  our  results 
for  the  impedance  of  the  SSM  and  ILSV  models  described  previously.  . The 
models  described  in  Belytschko,  et  al.,  (1976)  did  not  match  the  experimental 
impedance  satisfactorily,  so  the  alterations  necessary  to  match  experimental 
impedance  are  described  in  these  sections.  In  addition,  these  impedance  studies 
were  used  to  identify  the  stiffnesses  and  damping  constants  of  the  elements  inter- 
connecting the  viscera  and  vertebrae . 

Til . 2 Mathematical  Preliminaries 

A transfer  function  for  a linear  system  is  the  ratio  of  the  transforms  of  an 
input  variable  and  an  output  variable.  Hence,  impedance  is  a transfer  function. 

The  systems  we  are  considering  are  stable,  that  is,  they  are  systems  whose  natural 
responses  decay  with  time.  The  excitations  we  are  considering  are  such  that  they 
can  be  represented  by  exponentials  whose  frequencies  lie  along  the  imaginary  axis. 
In  this  case,  the  impedance  is  defined  by 

(3.1) 


Z(w) 


V(w) 


where  an  uppercase  letter  indicates  the  Fourier  transform  defined  by 


I 
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tf 
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f (t)  dt 

‘’o 

£(0 

J F(»)  (P® 
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(3.2) 
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Aa  an  example  consider  a simple  damped  oscillator  excited  by  a unit  step 
velocity -at  its" support-  The  equation  of  motion  fox  this  system  is 


nu^  + 0*2  + kx2  **  c^j.  + » 


which,  since  = 1 and  x^  • t can  be  written  as 


n#2  + c*2  + ^2  = c + kt 


The  solution  to  (3.5)  is 
-u&t 


■5-  Bln  8dc  + 1 


where 


^ *»  damping  ratio 


2 Jta 


(3.4) 


(3.5) 


(3.6) 


3 “ natural  circular  frequency  « JkJm 
p^  = damped  natural  circular  frequency 

= 0 Jl  - IA2 

The  force  transmitted  to  the  support  is  given  by 
f(t)  = k(x2  - x^  + c(^2  - x^)  = -mx2 

or,  using  the  solution  for  the  step  velocity  input  given  by  Eq.  (3.6), 

r 2 - m . 1 


(3.7) 


|(l-2»*2)  sin  P^t  + 2p  cos  p^tj-  . (3.8) 


The  transforms  required  in  Eqs.  (3.1)  are 


.03 

F(ib)  = J £(t)  e dt 


(3.9) 


(jiii)2  + 2^PJw  + g2 


• ‘ .Z<- 


i 


' . •»  ---  •-  i?r3,  • fv  -c,... 

j i*f>  £ * H'.t:  2?yt3s\'t  v^w'-  ( t2  fi  Z9- 


■r:r,:'i  ; , . £.  i 1 i$  ■ •>  .'  *'  'f  - ^ - 


V(u»)  »=  | dt  • " 


' : (3.10)  ; 


. J»  M * . 


The  impedance  of  this  system  is  then  given  by  the  magnitude  of  the  ratio  of 


Eqs.  (3.9)  and  (3.10): 


m &i)J  &2  4-  (2ptu)Z 


v/o2  - w2)  + <2 


(3. XI) 


If  we  define  the  frequency  ratio  r = fc/g,  Fq.  (3.11)  can  be  written  ae 


l/i+^pr)' 


7(l.-rZ)  + (2pr)2 


(3.12) 


For  negligible  damping  (p  s«Q),  2(r)  -•  « at  r “■  1,  or  » » g,  that  is,  the 


impedance  becomes  infinite  for  the  undamped  system  when  the  driving  frequency  is 


equal  to  the  natural  frequency.  For  increasingly  larger  p,  the  maximum  of  (3.12) 


occurs,  further  and  further  to  the  right  of  r » .1  and  the  maximum  impedance 


decreases. 


As  the  number  of  degrees  of  freedom,  and  heaee  complexity,  of  a mechanical 


system  increases,  it  becomes  increasingly  difficult  to  obtain  analytical 


expressions  for  the  response  sad  hence  the  impedance  of  the  system.  A time  history 


response  obtained  by  a transient  analysis  of  the  system  is  therefore  used  to  obtain 


the  impedance.  The  transient  raspoase  is  in  the  form  of  discrete  values  at  equal 


time  intervals.  If  tha  multidggree  of  freedom  system  being  considered  ia  linear 


and  stable,  and  if  che  excitation  is  expressible  in  tarsus  ai  exponential  components 


4m  £ 

of  the  form  eJ  n , where  u>n  now  represents  a discrete  spectrum  of  frequencies,  the 
response  can  be  transformed  by  means  of  the  discrete  Fourier  transform  (DFT)  vie 
its  computer  based  algorithm,  *"v>e  Fast  Fourier  Trans £orm  (FFT)  (Brigham  (1974) , 
Holmes’ (1976)). 

The  DFT  form  of  the  Fourier  transform  pair  inKqs.  (3.2)  and  (3.3)  is 
N-l 

F(wn)  - £ f(kfit)  e'H  Sc  (3.13) 

k=0 
N-l 

f(kAt)  = F(u»n)  ejwn  Sc  ’ (3.14) 

n=0 

where  At  is  the  sampling  interval,  i.e.  the  time  step  used  in  the  transient 


analysis,  N the  number  of  steps,  and  u>n  the  nth  frequency,  which  is  given  by 


_ 2rn 


(3.15) 


The  impedance  of  a model  is  obtained  by  prescribing  a unit  step  velocity  at 
the  driving  point,  which  in  the  simulations  reported  here  is  either  the  seat  or  the 
pelvis.  The  force  time  history  at  the  driving  point,  f(t),  and  its  discrete  Fourier 
transform  are  then  obtained.  The  transform  of  the  velocity  time  history  is  given 
by  Eq.  (3.10),  and  the  ratio  of  the  magnitudes  of  the  two  transforms,  Eq.  (3.1),  gives 
the  driving  point  impedance  that  is  used  here;  phase  angles  were  not  determined. 

A drawback  of  this  approach  is  that  a numerical  solution  of  a system  with  a 
unit  step  velocity  applied  directly  to  a mass  (e.g,  the  pelvis)  does  not  take  into 
account  the  mass  at  which  the  velocity  is  applied  because  tbs  acceleration  of  the 
mass,  which  is  a unit  impulse  analytically,  is  sero  for  the  discrete  system  except 
ot  the  time  aero,  where  it  is  not  computed.  In  other  words,  the  acceleration  time 
history  associated  with  a unit  step  velocity  is  an  impulse  function,  and  this 
impulse  will  not  appear  in  the  numerical  results.  Hence  the  internal  force  at  the 
mass  will  not  reflect  the  presence  of  the  mass. 
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However,  it  should  be  noted  that  when  the  system  is  driven  through  the  mass, 

the  mass  is  essentially  in  parallel  to  the  rest  of  the  system.  It  is  therefore 
possible  to  compute  the  impedance  of  the  entire  system  by  obtaining  .the  impedance 
of  the  system  without  the  mass  as  described  above,  and  then  adding  to  thistthe 
impedance  of  the  mass.  / ’ *-; 

This  difficulty  occurs  only  when  the  driving  point  is  the  pelvic  massy  However 
when  the  seat  is  the  driving  point,  the  seat-buttocks  stiffness  dominates  the 
impedance.  Since  it  is  experimentally  almost  impossible  to  drive  a human,  being  at 
the  pelvis,  the  impedance  curves  which  do  not  include  the  seat  are  physically 
unrealizable,  but  they  are  indicative  of  the  frequency  character  of  the  spine  itself 
In  Appendix  A,  the  DFT  method  used  here  and  the  programming  details  are 
described. 

III. 3 Suinnary  of  Impedance  Literature 

In  this  section  we  will  review  the  literature  on  experimentally  determined 
impedances  of  human  subjects.  The  papers  oj:  msjor  interest  are  those  of  Coermsnn, 
et  el.,  (1959-1960)  and  Vogt,  et  al,,  (1968).  Coermann  determined  the  impedance  of 
sitting  and  standing  subjects  by  shake  table  induced  longitudinal  vibrations  through 
a frequency  range  of  1 to  20  Hs.  Figure  3.1  shows  the  median  impedance  Of  five 
sitting  subjects  reported  by  Coeraann.  Also  shown  is  an  impedance  curve  for  a one 
degree  of  freedom  system,  Eq.  3.12,  with  Coermann’s  values  for  the  mass,  stiffness, 
and  fraction  of  critical  damping.  Figure  3.2  shows  the  effects  of  variation  in  body 
postures  as  determined  experimentally  by  Coarmann.  The  impedance  of  the  abdominal 
viscera  was  also  studied  by  Coertuaan  by  placing  the  subject  in  a supine  position  on 
a shake  table  and  subjecting  him  to  longitudinal  harmonic  vibrations  through  a 
frequency  range  of  1 to  15  Hs.  Motion  of  the  skeleton  was  suppressed  by  using  metal 
brackets  at  the  head,  shoulders  and  feet,  to  clamp  the  subject  to  the  table. 
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Figure  3.1.  Comparison  of  the  median  of  the  experimentally  measured  impedances  of 
5 sitting  subjects  (longitudinal  vibrations)  to  a single  degree-of- 
freedom  system  with  m ■ 84000  grams,  0 * 10tt  and  ^ •>  0.312.  (Coeraann, 
et  al.,  (I960));  (1)  experimental;  (2)  computed. 
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Figure  3.2,  Influence  of  body  posture  on  the  mechanical  impedance  of  a human  subject; 
(1)  sitting  relaxed;  (2)  sitting  erect;  (3)  standing  erect. 

(Coertaaim,  et  ai. , (I960)). 
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Figure  3.3  shows  the  variation  of  abdominal  wall  displacement  with  excitation 
frequency  as  measured  by  Coeroann,  et  al.  A resonance  in  the  abdominal  wall  dis- 
placement was  detected  between  3 and  4 Hz,  so  this  would  indicate  an  impedance  peak 
at  between  4 and  5 Hz.  Restricting  the  mobility  of  the  abdomer  resulted  in  a shift 
of  the  abdominal  wall  displacement  peak,.tp  between, 6 and  7 Hz. — - ■ 

Vogt,  et  al.,  (1968)  measured  the  impedance  of  ten  male  subjects.  The  -subjects 

S • • ‘ 

were  seated  and  loosely  restrained.  Figure  3.4  shows  the  f**era ge  impedance  of 
these  subjects.  An  important  characteristic  of  the  above  experimental  irapa^arcee 
is  the  peak  which  occurs  around  5 Hz.  A second  peak,  which  is  generally  considerably 
smaller  than  the  first,  occurs  around  14  Hz.  These  peaks  in  the  impedance  curve 
indicate  the  excitation  frequencies  at  which  the  maximum  energy  is  transferred  to 
the  body.  Coermann  suggests  that  the  5 Hz  peak  is  caused  by  resonant  mot ion  Of  the 
upper  torso  in  connection  with  the  bending  elasticity  of  the  pelvis  and  spine*  and 
that  the  14  Hz  peak  is  probably  attributable  to  another  elasticity  of  the  pelvis. 

Vykukal  (1968),  exposed  4 subjects  in  a semi-supine  position  to  a vertical 
acceleration  of  ± 0.4  G in  a frequency  range  from  2^  t-  20  Hz,  combined  with 
linear  acceleration  of  1,  2^,  and  4 6.  The  mechanical  impedance  of  each  subject  was 
recorded.  He  observed  that  for  the  higher  linear  accelerations,  the  stiffness 
increased,  the  damping  decreased,  and  the  overall  impedance  increased.  The 
resonances  at  higher  frequencies  became  more  predominant  in  magnitude  because  of 
increased  coupling  of  the  boc’y  systems  when  subjected  co  a high  G environment. 
Unfortunately,  neither  the  direction  of  the  harmonic  oscillation  vector  nor  the 
direction  of  the  linear  ac  :eleraciou  vector  with  respect  to  the  subject's  longitudinal 
axis  were  specified. 

Woods  (1967),  reported  transmissibilities  for  longitudinal  and  lateral  vibrations 
of  three  subjects  in  the  frequency  range  of  1 to  10  Hz.  Results  are  also  reported 
for  the  effect  of  random  vertical  and  lateral  vibrations.  All  of  the  transmissibility 
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Figure  3.3.  Abdominal  wall  displacement  versus  excitation  frequency  for  a human 
subject  in  the  supine  position  exposed  to  longitudinal  vibrations. 
(Coerva&nn,  ot  al,,  (i960)),  . 
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Comparison  of  experimentally,  ane lyrically  and  numerically  determined 
impedance  curves;  (1)  experimentally  determined  impedance  curve  of 
Vogt,  et  el.,  (2)  analytically  computed  impedance  curve  for  single 
degree  of  freedom  system;  (3)  numerically  determined  impedance  curve 
for  single  degree  of  freedom  system. 
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curve'  for  the  sinusoidal,  longitudinal  oscillations  had  a peak  between  4 to  6 Hz. 

A peak  in  the  trensmissibiiity  for  the  sinusoidal,  lateral  oscillations  occured 

t 

at  1.5  Hz.  The  vertical  transmisflibility  curves  ye re  obtained  by  plotting  the  ratio 
of  the  accelerations  of  the  shoulder  to  the  acceleration  of  the  seat  and  the  accel- 
eration  of  the  thigh  to  the  acceleration  of  the  seat  as  a function  of  the  input 
frequency.  The  lateral  transmissibility  curves  were  obtained  by  plotting  the  ratios 
of  hip  acceleration  to  the  seat  acceleration,  the  knee  acceleration  to  the  seat 
acceleration,  the  shoulder  acceleration  to  the  seat  acceleration,  and  the  head 
acceleration  to  the  seat  acceleration  as  a function  of  the  input  frequency.  All  of 
these  curves  had  a peak  at  1.5  Hz.  It  is  not  clear  whether  these  peaks  are  a con- 
sequence of  human  body  response  or  the  subjects'  interaction  with  the  seat  and 
restraint  system.  However,  the  reported  peak  in  the  4 to  6 Hs  range  correlates  well 
with  the  results  reported  in  Coermann,  et  al.  and  Vogt,  et  al. 

It  is  also  worthwhile  to  consider  the  impedance  characteristics  of  previous 
models.  Payne  (1972)  developed  a four  degree  of  freedom  model  of  a seated  man. 

It  consisted  of  rigid  masses  representing  the  head,  upper  torso,  abdominal  viscera 
and  pelvis,  interconnected  by  springs  and  dampers.  A spring  and  a damper  were 
included  to  represent  the  elastic  and  viscous  properties  at  the  buttocks.  Impedance 
curves  for  this  model  exhibit  a peak  at  9 Hz  and  a severe  dip  in  the  impedance  curve 
at  12  Hz  but  no  second  resonance  point,  so  they  do  not  agree  well  with  Vogt,  et  al. 
or  Coermann,  et  al. 

The  impedance  of  the  model  in  AHRL-TR-76-10  was  not  given,  but  modal  analyses 
were  performed.  These  modal  analyses  showed  that  the  ligsmented  spine,  with  head 
and  pelvis  but  no  separate  elements  for  the  viscera  and  buttocks,  aad  with  only  axial 
motion  allowed,  yielded  a lowest  natural  frequency  of  17  Hz.  The  mode  shape 
associated  with  this  frequency  is  approximately  that  of  a free- free  rod.  This  result 
Indicates  that  tLa  5 Hz  peak  in  the  impedance  curve  reported  experimentally 
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represents  a resonance  of  either  the  viscera , the  buttocks*  or  the  beading  response 
of  the  spine.  The  relatively  low  mass  of  the  viscera  would  tend  to  rule,  it  out  as  , ,! 
the  sole  source  of  this  peak.  . It  will  be  shown  subsequently  that  the  first  peak 
in  the  impedance  curve  depends  strongly  on  the  stiffness  of  the  buttocks  and  seat.  . 

It  will  also  be  shown  that  if  the  possibility  of  parametric  excitation  of  the 
flexural  spine  modes  is  considered,  a 5 to  7 Hz  peak  is  present  in  the  impedance 
curve  even  without  the  representation  of  the  buttocks,  but  its  magnitude  la  smaller 
than  that  found  experimentally.  Thus,  the  first  resonance  peak  observed  experimen- 
tally appears  to  arise  from  several  components  of  the  human  body:  the  viscera, 
the  seat-buttocks,  and  the  bending  response  of  the  spine. 

III. 4 Impedance  Results 

Impedance  curves  were  obtained  for  two  of  the  models,  the  simplified  spine 
model,  SSM,  and  the  Isolated  ligamentous  spine  with  viscera,  ILSV.  These  impedance 
studies  were  used  to  obtain  the  stiffnesses  and  damping  constants  for  the  elements 
interconnecting  the  spine  and  viscera.  The  stiffnesses  were  first  determined  in 
the  SSM  model  by  parametric  studies. 

The  damping  of  these  elements  was  then  chosen  so  that  the  lowest  "lateral 
visceral  mode"  (i.e.  motion  of  the  viscera  laterally  relative  to  the  spine)  was 
damped  by  the  same  amount  as  Che  lowest  axial  visceral  mode.  Based  on  the  stiff- 
nesses of  the  SSM  viscera/spine  interconnect  lag  elements,  initial  values  for  the 
ILSV  interconnecting  elements  were  chosen.  Only  minor  adjustments  were  necessary 
to  match  the  experimental  impedance. 

Prior  to  applying  the  FFT  impedance  determination  method  to  the  spina  models, 
it  was  applied  to  a one  degree  of  freedom  model  using  the  mass,  damping  and  frequency 
values  used  by  Coemaan,  et  al.  (I960)  for  his  one  degree  o£  freedom  model  of  a 
185  lb  man.  The  numerical  results  are  shown  in  Fig.  3.4  along  with  Che  analytical 
results  of  Eq.  (3.11).  The  duration  of  the  numerical  simulation  of  the  step 
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response  hers  was  1.02  anc,  which  involved  256  tin*  seeps  of  4 milliseconds.  As  can 
be  seen,  the  results  compare  very  well. 

We  will  now  give  some  of  the  results  for  the  SSM  model-  The  parametric  studies 
used  to  determine  the  spine-viscera  connection  elements  are  not  reported;  the  major 
emphasis  in  the  following  studies  is  the  origin  of  the  various  resonances  and  the 
influence  of  certain  parameters,  such  as  seat-buttocks  stiffness,  on  the  impedance. 

Figure  3.5  compares  the  impedance  curves  for  the  SSM  without  buttocks  to  the 

averaged  experimentally  obtained  curve  of  Vogt,  et  al.  The  latter  were  obtained  by 

vertical  shaking  of  loosely  restrained  seated  subjects,  so  the  buttocks  must  be 

evidenced  in  the  observed  impedance.  The  pelvic  mass  here  is  1.62  x 10  grams 

4 

(36  lbs),  but  the  total  mass  of  5.56  x 10  grams  (123  lbs)  does  not  include  the 

4 

upper  or  lover  leg  mass.  The  average  mass  of  Vogt's  subjects  was  6.94  x 10  grams 
(153  lbs)  (without  the  lower  legs  and  feet,  which  ore  not  included  because  the 
footrest  of  the  experiment  was  not  harmonically  excited). 

The  impedance  of  the  spine  without  buttocks,  as  shown  in  Fig.  35,  exhibits 
peaks  at  approximately  5.70  Hz  and  12.70  Hz  respectively,  both  with  and  without  the 
pelvic  mass.  In  order  to  more  clearly  determine  vhiwh  subsystems  are  causing 
these  peaks,  the  impedances  of  the  spine  and  the  viscera  are  given  in  Fig.  3.6.  The 
axial  response  of  the  spine  alone  has  one  peak  at  10.50  Hz,  while  the  visceral 
column  haa  a peak  at  approximately  5.37  Hz.  These  agree  with  the  lowest  frequencies, 
obtained  from  an  axial  modal  analysis  of  the  spinal  column  and  visceral  column  with 
both  ends  fixed.  The  impedance  of  the  spine  and  viscera  combined,  though  dominated 
by  the  peak  at  approximately  10.74  Hs  (again  the  first  fixed-free  spine  axial  mode) 
also  exhibits  the  presence  of  the  viscera  with  the  "beginnings  of  a peak"  at 
approx leas  tel y 5,70  Hz.  this  curve  is  dominated  by  the  spine  peak  because  there  is 
considerably  more  ’nass  associated  with  cite  spine  (e.g  head,  arms,  ete.\  than  with 
the  viscera.  Nevertheless,  the  visceral  contribution  to  the  impedance  curve  is 
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Impedance  curve*  tc  the-  SSM  without  buttocks;  (1)  experimental  curve 
of  Vogt,  et  el.,;  (2)  S.‘-  i pelvic  mass  impedance ; (3)  SSM  total  sass 
impedance;  (4>  isped^nct  of  SSM  without  pelvic  oase;  (S)  impedance 
of  SSM  with  pelvic  * se. 
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It  is  of  interest  to  note  that  upon  dropping  the  restriction  of  .axial  motion, 
we  obtain  the  curves  Fig.  3.5,  which  have  a peak  at  12.70  Hz.  In  the  modal  analysis 
results  of  SSM,  the  lowest  visceral  mode  and  the  third  spine  bending  mode  (in  which 
most  of  the  energy  is  transmitted  to  the  lumbar  region)  were  at  5.25  H2,  -The 
discrete  internal  force  time  history  results  for  the  spine  element  between  Si  and  O 
L3  were  in  fact  dominated  by  this  frequency.  This  indicates  that  the  third  ...... 

bending  mode  has  been  parametrically  excited.  A variation  %f  the  bending  damping 
or  visceral  damping  had., its  greatest  ef feat  JUl  the  region  of  ,the  5.70  Hz  peak.  v ?•’. 
These  results  indicate  that  SSM  has  a.  5.70  Hz  impedance  peak  due  to  the  first 
visceral  mode,  and  parametric  excitation  of  the  third  bending  mode,  even  in  the 
absence  of  a seat-buttocks  representation. 

In  determining  the  effects  of  the  buttocks,  the  buttocks  were  modeled  by  a 

spring  and  damper  between  the  pelvis  centroid  and  the  seat.  Payne  ;('1972)  gives  a 

value  of  6.55  x 10^  dynea/cm  (374  Ib/in)  for  the  buttocks'  stiffness.  Buttocks 

7 8 

stiffnesses  of.  6.55  x 10  stnd  1.31  x 10  dynes/cm  were  used  here  and  half  of  the  up - 

4 

per  leg  mass  was  added  to  the  pelvic  mass  t.o  yield  a total  SSM  mass  of  6.88  x 10 
grams  (152  lbs),  because  in  the  experiment  the  upper  leg  mass  was  vibrated  along 
with  the  upper  body. 

Figure  3.7  .-.hows  the  impedance  of  SSM  with  these  two  values  of  buttocks 
stiffness.  Both  curves  exhibit  a peak  at  4.30  Hz  and  a much  smaller  peak  at  9.30  Hz, 
with  a slight  shift  to  the  right  for.  the  stiffen  buttocks  model.  Modal  analysis 
shows  that  the  first  axial  mode  of  the  buttocks-spine  series  combination  occurs  at 
4.10  Hz  and  consists  essentially  of  rigid  body  motion  of  the  spine  relative  to  the 
seat.  The  9.30  Hz  peak  appears  to  correspond  to  a 9.10  Hz  mode  which  consists 
primarily  of  the  second  visceral  mode, 
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Figure  3.7.  Impedance  of  SSM  with  buttocks:  (1)  experimental  curve  of  Vogt,  et  al., 

(2)  and  (3)  SSM  with  buttocks  stiffnesses  of  6.55  x 107  and 
8 

1.31  x 10  dyne/cta  respectively. 


From  these  results  it  becomes  apparent  that  an  impedance  curve  obtained  by 
driving  the  buttocks  is  dominated  by  the  steatopygic  mode.  This  agrees  with  the 
findings  of  Payne  and  Band,  (1969),  who  comment  that  the  variability  in  impedance 
measurements  of  seated  human  subjects  is  very  likely  due  to  variations  In  buttock 
size,  muscle  tension,  and  tonus.  It  follows  from  this  that  impedance  measurements 
of  this  type  do  not  result  in  accurate  information  as  to  resonances  of  the  sub- 
systems.  However,  it  is  interesting  to 'observe  that  the\ shift  in  impedance  peaks 
due  to  a twcfold  change  in  buttocks  stiffness  is  no t\ severe. 

Although  not  apparent  in  bis  results,  Vogt,  et  -al.-,  (1968),  reported  that  in 
most  cases,  a peak  around  10  Hz  is^.also  evident.  As  discussed  above,  the  impedance 
curve  for  the  SSM  with  the  buttocks  exhibits  a small  peak  around  10  Hz.  However  it 
cannot  be  established  whether  the  peak  in  Vogt,  et  al  is  caused  by  the  same  mechanism 

\ } i 

which  causes  the  10  Hz  peak  in  the  SSM  curves.  Their  impedance  a Iso  exhibits  a email 
peak  at  14  Hz.  This  peak  does  not  show  up  in  the  SSM  curves. 

It  is  clear  from  these  results  that  the  SSM  model  duplicates  experimentally 
determined  impedances  very  well.  It  not  only  replicates  the  peak  magnitudes  and 
the  locations  of  these  peaks,  but  it  gives  considerable  insight  into  the  resonance 
mechanisms. 

We  will  now  consider  the  impedance  studies  of  the  xiore  complex  model,  ILSV, 
isolated  ligamentous  spine  with  viscera,  which  was  described  in  Section  II. 2. 

Figure  3.8  shows  the  impedance  curves  for  the  ILSV  without  buttocks.  For  comparison, 

the  experimental  curve  of  Vogt,  et  al.,  (1968)  is  also  shown.  The  pelvic  mass, as 

4 4 

for  the  SSM,  is  1.62  s 10  gratss  <36  lbs),  the  total  b&ss  is  4.98  x 10  grams 

(110  lbs).  The  total  mass  does  not  include  the  upphr  os  lower  leg  mass,  nor  the 

ana  mass,  whief  w»rc  included  in  the  3SK.  iapadancas  of  tha  ILSV  without  and 

with  the  pelvic  are 
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Figure  3.8.  Impedance  curves  for  the  isolated  ligamentous  spine  with  viscera 
(ILSV)  without  buttocks:  (l)  expevime  al  curve  of  Vogt,  et  el.; 

(2)  ILSV  pelvic  mass  impedance;  (3)  ilsv  total  maes  impedance; 

(4)  ILSV  impedance  without  pelvic  mass;  (51  ILSV  impedance  including 
pelvic  mass. 
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The  impedance  exhibits  peaks  at  approximately  6.0  and  13.5  Hz,  which  agrees 
well  with  the  corresponding  SSM  peaks  at  approximately  5.7  and  12,7  Hz.  The  13,5 
Hz  peak  corresponds  to  the  first  spine  axial  mode,  13.17  Hz,  obtained  in  an  ILSV 
modal  analysis  with  the  pelvis  fixed.  The  6.0  Hz  peak  results  from  a combination 
of  the  first  visceral  mode  and  what  appears  to  be  a bending  mode.  The  visceral 
frequency  is  5.5  Hz,  and  there  are  three  bending  modes  in  the  5-6  Hz  region. 


The  impedance  results  for  the  ILSV  with  buttocks  are  shown  in  Fig.  3*9.  The 


same  buttocks  stiffness  as  in  the  SSM,  6.55  x 10/  dynesVcm,  (374  Ib/in)  wan  used. 


Half  of  the  upper  leg  mass  was  added  to  the  pelvic  mass  to  yield  a total  ILSV  mass^" 


of  6.26  x 10  grams  (138  lbs).  The  impedance  of  ILSV  is  given  with  buttock  damping 
ratios  of  0.3  and  0.15,  res  pec  Lively.  TheA.8  jgz  peak  corresponds  to  the  steatopygic 
mode  (4.4  Hz)  and'as  in  the  SSM,  this  mode  dominates  the  impedance  curve.  A small 
peak  is  also  observed  at  approximately  10.0  Hz,  which  appears  to  be  associated  with 
the  second  visceral  mode. 

Figures  3.10  and  3.11  compare  the  SSM  and  ILSV  impedance  curves  without  and 
with  buttocks,  respectively.  Of  primary  interest  is  the  similarity  of  the  SSM  and^;. 
ILSV  curves.  This  indicates  that  SSM  and  ILSV  would  respond  almost  identically 
for  dynamic  input  with  a frequency  content  of  16  Hz  or  less.  . 
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Figure  3.9. 
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| Figure  3,10.  Comparison  of  impedance  curves  for  SSM  and  XLSV  without  buttocks 

1 (1)  experimental  curve  of  Vogt,  et  al.;  (2)  and  C3)  SSM  and  ILSV 

\ without  pelvic  masses  respectively;  (4)  and  (5)  SSM  and  1LSV 

f with  pelvic  masses  respectively. 
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Figure  3.11.  Comparison  of  impedance  curves  for  ecu  n ^ „„„ 

(1)  experimental  curve  of  Vogt,  et  al  mT  ^ buttocks; 

<3>  l«p.da„ceB„f  ILSV  with  bitwck.!"6  °f  SSM 
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A simulation with  the  spine  head  model  yieldf  as  part  of  the  results. 


the  normal  force  N,  and  bending  moment,  M,  time,  his  tor  tea  at  each  vertebral 


level.  In  this  Chapter,  the  evaluation  of \ injury  [potential  as  a function  . 


of  N and  M will  be  considered.  We  will  call  this  ! function! the -Injury  Potential 

' ‘ „.,.J  / \ T 

Function.  1 


Mathematically,  we  can  express  our  task  as  the  determination  of  where 


IPFj  - £<  lMj  1 ■ bjl*  Alj-42^ 


(4.1) 


d ; ■ ■ 


in  which  IPF^  * injury  potential  function 


. i 


|m,},  Jn . J =*  bending  moment  and  normal  force  magnitudes  when 

J r ' \ 


IPp  is  maximum  in  vertebral  body  j 

Y 


Alj,  A2^  * parameters  representing,  the  failure  criteria  of 

vertebral  body  j. 

Belytschko,  et  al.,  (1976)  discussed  an  approach  for  the  preliminary 
evaluation  of  injury  potential.  They  idealized  the  vertebral  body  as  a 
cylindrical  shell  of  radius  rQ  and  height  h (see  Fig.  4.1).  The  shell 
consisted  of  cortical  or  compact  bone  with  an  interior  (vertebral  core)  of 
trabecular  or  soft  bone  of  radius  r^. 

For  pure  axial  compression,  the  stresses  in  the  cortical  shell  and 


vertebral  cere  are  related  by 

sr,  - (4.2) 

and  equilibrium  requires  that 


<7.  “ (R  “ <?  A)  /A 


(4.3) 


where 


°i,0o  = stress  in  trabeeul^r  and  cortical  botte.  reapecttvely 

Ei»E0  3 Young's  modulus  in  trabecular  and  cortical  bone,  respectively 
N = applied  axial  force 


Aq  = area  of  cortical  shell  » n(r2  - r 2) 

o X 

2 

“ area  of  vertebral  core  = rir^  • 


Eqs.  (4.2)  and  (4.3)  yield  ao  expression  for  the  stress  in  the  cortical  shell. 


"fcirt2  + Vr=2  ri2)1 


(4.4) 


From  a compression  test,  we  can  obtain  an  axial  force,  N .which  is  the 

axial  force  that  can  be  sustained  by  a vertebra.  The  corresponding  stress, 

as  given  by  Eq.  (4.4),  is  then  the  maximum  safe  stress  for  the  cortical  bone. 

hQS  been  determined,  we  can  use  it  to  obtain  H the  moment 

max, 

which  would  produce  a stress  in  the  cortical  shell.  The  flexure 

formula  from  simple  beam  theory  results  in  the  relationship 


4Mr  E 
o o 


r 4 4 A * 


(4.5) 


from  which 


M • 
max 


"[E.r,4  + S U 4 

11  O 0 

4r  E 

O O 


- «>o  - « 
‘i  )J  *o 


(4.6) 


For  combined  bending  and  axial  loading,  <r  ***  is  given  by  the  superposition 


of  Eqs,  (4.4)  aud  (4.5), 


a + s . i 

m * w * # 

max  max 


(4.8) 


The  above  formulation  lends  itself  directly  to  an  injury  potential  criterion, 


IPF 


I", I + hi  <1 


Al. 


A2 . 


(4.9) 


j j 

where  Jm^ | and  jN^ | are  as  defined  in  Eq.  (4,1),  and 
j ® vertebral  level 

Al^  ® maximum  bending  moment  in  pure  bending 

A2j  9 maximum  compression  force  in  pure  compression. 

Eq.  (4.9)  was  developed  by  considering  the  stresses  in  the  vertebral  body's 
cortical  shell.  No  direct  attention  was  given  to  the  stresses  in  the  vertebral 
core.  The  elastic  modulus  of  the  vertebral  core  is  several  orders  of  mafpiitude 
lower  than  that  of  the  cortical  shell.  Hence,  the  stresses  will  be  several 
orders  of  magnitude  larger  in  the  cortical  shell  even  under  pure  axial  loading. 
When  we  consider  bending  stresses  also,  this  disparity  becomes  even  greater. 

It  is  this  large  difference  in  cortical  shell  and  vertebral  core  stresses, 
particularly  under  the  combined  action  of  bending  and  axial  loading,  which 
justifies  the  use  of  the  cortical  shell  stress  as  the  basis  for  the  injury 
potential  criteria. 

Eq.  (4.9)  is  essentially  the  same  relationship  used  by  Belytschko,  et  al., 
(1976)  in  their  preliminary  evaluation  of  injury  potential.  However,  rather 
than  determining  the  maximum  of  IPF^  for  the  entire  time  history,  they  used  the 
individual  maxiumus  of  {n^|  and  |n ^ j to  determine  IPF^.  This  leads  to  larger 
values  of  IPF^- 

Before  the  injury  potential  function,  Eq.  (4.9),  can  be  put  to  use.  It  is 


necessary  to  obtain  values  for  the  vertebral  material  properties  and  geometries 


IV. 2 Compressive  Strength  of  Vertebrae 

Payne  (1972)  examined  the  work  of  Geertu  (1946)  in  which  Che  breaking  loads 

of  vertebra*8  T8  through  L5  wt re  determined  £os  subjects  varying  ir  age  from 

19  to  46.  These  are  reproduced  in  Table  4.1.  Also  included  in  Table  4.1 

are  three  values  of  breaking  atrergths  (T8,  T12,  L2,  Ll)  determined  by  Crocker 

and  Higgins  (1966).  The  average  lumbar  vertebral  breaking  load  of  the  vertebrae 

listed  in  this  table  is  9.28  x 1G®  dynes  (2090  lbs.).  Brown,  et  &1,»  (1957) 

reported  an  ultimate  compressive  load  varying  from  4.45  x 10  to  5.80  x 10° 

dynes  (1000  to  1300  lbs.)  for  5 lumbar  vertebrae  taker*  from  three  measurements. 

8 

Perey  (1957)  reports  an  average  value  of  5.88  x 10  dynes  (1320  lbs.)  for  the 

lumbar  vertebral  breaking  load.  Payne,  however,  comments  that  Perey's  failure 

to  "true-up"  his  vertebrae  to  insure  for  uniform  transmission  of  the  load  from 

the  testing  apparatus  to  the  specimen  resulted  in  induced  stresses  higher  than 

for  pure  uniform  compression,  thus  resulting  in  e lower  apparent  breaking  load. 

Kockoff,  et  al.,  (1969)  determined  the  maximum  non-destructive  compressive 

strengths  (elastic  limit)  of  50  lumbar  vertebrae  from  the  spines  of  32  fresh 

cadavers.  The  subjects  included  both  males  and  females  and  spanned  an 

approximate  age  range  of  2d  to  80  years.  Rockoff  reports  chat  all  of  the 

specimens  from  subjects  of  age  less  than  40  years  had  strengths  greater  than 
7 2 

800  ps<  (5.52  x 10  dyna/cm  ) while  essentially  all  of  the  specimens  greater 
than  40  years  had  strengths  less  than  800  psi.  If  we  assume  this  is  the  average 
strength  of  the  specimens  and  assume  a lumbar  vertebral  body  cross-sectional 
area  of  16-51  cn~  (from  Payne) , ve  obtain  an  average  elastic  limit  load  of 

Q 

9.11  x 10  dynes  (2050  lbs).  Rockoff  also  considers  the  rel*  v . ‘it 
limit  strength  contribution  oi  the  vertebral  cortex  and  the  trabecular  bone, 
and  the  variation  of  the  elastic  limit  with  the  ash  content  or  ihysical  density 
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Table  4.1,  Breaking  loads  of  vertebrae  T8  through  L5 
a a determined  by  Geertz,  and  Crocker  and 
Higgins,  reported  by  Payne 


Age 

[(Years) 

[Geertz 


Breaking  Load  in  Dynes  x I0£ 


19  

_21 6.3  7.8 

21  

J_3 5.3  6.0  

33 6j3_ 

J6 5.9  7.1  7.5 

38 


TlO  Til 

T12 

1 "** 

Ll  L2  L3  U 

L5 

7.4 

'! 

t 

7,1  8.8 

7-8  | 

1 8.8 

L 9.7  , 

|10.0| 

7.1 

! 6.8 

r 8.2I 

i 

i 

7.8  i 

7.8 


10.8 

S.g- 

9.3 


Crocker 

and 

Higgins 

47 5.3 

m.Jxi'-1'*  Is, 8 


9,5  I 11.0  | 

6,7  ! 7,3  7.6  7. si  7,9  8.5  ! 9,6  9.6  l 

sjjsb^r  vertebrae  9.2 


i , ^ \ 

iK-fi-sr, ' for  lumbar  vertabree 


data  on t y 


Ol;'A  w 


of  the  trabeculate  bone.  He  found  'chat  the  cortex  generally  contributes  45 
to  75%  of  the  vertebral  body  strength  regardless  of  the  ash  content  of 
the  trabecular  bone. 

These  sources  as  well  r.s  several  others  indicate  that  the  upper  limit  of 

9 

lumbar  vertebral  breaking  strength  is  approximately  1.0  x 10  dynes  (2250  lbs.). 
Although  Rockoff  reported  that  his  strengths  were  elastic  limit  loads,  it  appears 
that  the  load  deflection  curves  for  vertebra,  are  essentially  linear  up  to  the 
point  of  failure.  Since  the  elastic  limit  is  defined  as  the  point  at  which 
the  load-deflection  curve  deviates  from  linearity,  it  is  quite  close  to  the 
breaking  strength. 

,.V , 3 Elastic  Moduli  of  the  Vertebral  Cortex  and  the  Trabecular  Bone 

Data  for  the  elastic  modulus  of  the  vertebral  cortex  could  not  be  found 
in  the  literature.  However,  several  investigators  have  reported  elastic 
moduli  of  cortical  bone  from  other  bones.  McElhaney  (1966)  reported  a value 

11  2 (j 

of  1.52  x 10  dyne/cm  (2.2  x 10  psi)  for  the  elastic  modulus  of  compact 

bone  from  an  embalmed  human  femur.  Evans  (1970)  reports  values  of  1,45  x 10 

1.59  x i0^~  and  1.74  x 10^  dyne/an“  (2.1  x 10^,  2.51  x 10^  and  2.52  x 10^  psi) 

for  the  elastic  moduli  of  compact  bone  from  a wet  embalmed  adult  femur,  tibia, 

and  fibula  respect iv^y.  Yamada  (1970)  obtained  a mean  value  of  1.04  x 10^ 

2 6 

dyne/cm  (1.51  x 10  psi)  for  the  elastic  modulus  of  wet  femoral  compact  bone 

8 2 

of  people  20  - 39  years  of  age.  He  also  reports  a value  of  8.82  x 10  dyna/cm 
(1.28  x 10  psi)  for  the  elastic  modulus  of  wet  trabecular  bone  from  human 
lumbar  vertebrae . 

IV, 4 Vertebral  Geometric  Properties 

Belytschko,  et  1.,  based  the  equivalent  radii  of  the  vertebral  bodies 
on  Lanier's  (1939)  mean  values  of  the  transverse  and  sagittal  diameters.  If  . 
we  assume  an  elliptical  vertebral  body  cross-section  and  the  transverse  and 
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sagittal  diameters  are  given  by  2a^  and  2b^  (vertebral  level  j)  then  the 

equivalent  radius  is  r.  B ,b  . , 

J oj  oj 

The  value  which  Belytschko,  et  al.,  used  for  the  vertebral  cortical  bone 
thickness,  0.3  ran,  was  based  on  Kulak  (1974).  Kulak  measured  the  cortical  bone 
thickness  of  sectioned  vertebrae  by  viewing  them  under  a microscope  through 
a measuring  eyepiece.  These  vertebrae  were  obtained  from  the  lumbar  region 
of  five  s-bjects.  The  mean  value  obtained  was  0.32  mm.  It  is  possible  that 
the  cortical  shell  thicknesses  of  vertebrae  in  the  thoracic  region  are  less 
than  0.3  mm,  but  since  these  values  have  not  been  reported  in  the  literature 
to  the  authors'  knowledge,  the  assumption  of  a constant  vertebral  cortical 
shell  thickness  of  0.3  mm  for  all  vertebral  levels  will  be  used. 

IV. 5 Failure  Criteria  for  the  Injury  Potential  Function 

Table  4.2  lists  the  values  of  experimental  axial  compression  failure  loads 
used  by  Belytschko,  et  al.  They  evidently  based  their  values  for  T8  through 
L-5  primarily  on  the  work  of  Geerts  (1946).  The  values  for  Tl  through  T7 

are  based  on  the  assumption  that  the  vertebral  cross-sectional  area; 

and  hence  breaking  strength,  decreases  linearly  with  vertebral  level.  This 

is  very  close  to  true  for  levels  T8  through  L5. 

For  the  elastic  moduli  of  the  vertebral  cortex  and  trabecular  bone,  they 

11  8 2 

cite  values  from  Evans  of  1.5  x 10  and  7.35  x 10  dyne7cra  respectively 
6 4 

(2.18  x 10  and  1.07  x 10  psi)  which  appear  to  be  in  reasonable  agreement 
with  thoae  found  by  other  investigators. 

Comparing  the  parameters (breaking  strength,  elastic  moduli,  etc.)  which 
Belytschko,  et  al.  used  in  their  evaluation  of  injury  potential,  with  those 
discussed  in  the  previous  section,  it  is  apparent  that  the  failure  criteria 


Table  4.2.  Maximum  axial  compression  for  pure  compression  (A2)  and  maximum 
bending  recent  for  pure  bending  (Al)  used  by  Belytscbko,  et  al. 
in  their  preliminary  evaluation  of  injury  potential 


Vertebral 

Level 


A2<3 4> 

Al<2> 

X 

1 

y'-N 

o 

o 

N-y* 

10.29 

10.19 

2.18 

9,80 

9.87 

2.03 

9.31 

9.18 

1.98 

8.82 

8.40 

1.95 

8.33 

7.72 

1.91 

7.84 

7.11 

1,84 

7.35 

6.37 

1.82 

6.86 

5.64 

1.81 

1 

6.37 

4.98 

1.78 

5.88 

4.41 

1.73 

5.58 

3.98 

1.74 

5.29 

3.53 

1.78 

4.90 

3.16 

1.71 

4.60 

2.84 

1.69 

4,21 

2,  .49 

1.62 

3.92 

2.24 

1.57 

3.62 

1.96 

1.54 

2. 


2.203 


2.154 


2.078 


2.019 


1.972 


1.883 


1.7181 


1 . 690 


1.619 


1.538 


1 .437 


1.388 


1.271 


1 .226 


1.162 


(1)  Equivalent  radius  (cm) 

8 

(2)  Maximum  bending  moment  ror  pure  compression  (dyne, /cm  x 10  ) 
computed  from  Eq.  4.6 

(3)  Experimentally  determined  compressive  breaking  load  (dyne  x 10^) 
from  Geertz,  reported  by  Payne 

2 9 

(4)  Maximum  vertebral  cortex  stress  (dyne/cm  x 10')  frost  Eq.  4.4 


(maximum  axial  compression  load  and  maximum  bending  moment)  computed  from  these 
parameters  is  acceptable  within  the  range  of  accuracy  permitted  by  the  available 
data.  Consequently,  the  present  injury  potential  study  incorporated  identical 
parameters  in  the  development  of  the  injury  potential  function  described  in 
Section  IV, 1. 

This  injury  potential  analysis  has  been  incorporated  as  part  of  the  head 

spine  simulation  program  as  follows.  During  each  step  of  the  time  integration 

process,  the  axial  force  N and  the  two  moments  M and  M are  obtained  for  each 

y z 

vertebral  level  j.  Since  the  program  treats  the  vertebrae  as  rigid  bodies  and 
the  discs  as  deformable  elements,  the  axial  force  and  moments  are  not  directly  avail- 
able, so  the  vertebral  forces  at  level  j are  considered  as  the  average  of  the  forces 
in  the  discs  above  and  below  the  vertebral  body;  i.e. 


My  “ W'Sy!  + V 

N * + V 

• «ixl  + A2K) 


(4.10) 


where  I and  K are  the  secondary  nodes  at  the  top  and  bottom  of  the  vertebral 
body  j.  Although  the  local  x axes  of  the  disc  abeve  and  below  the  vertebral  body 
may  differ  slightly  in  orientation,  this  nas  been  neglected.  The  local  y and 
z directions  in  each  disc  correspond  to  the  moments  in  the  frontal  and 
sagittal  planes,  respectively . 

The  injury  potential  function  is  then  evaluated  for  both  moments 

I"!  K! 


(4.11) 


IPP  « — ■+■  -X 

irl  Aij  A2^ 

i»l  \\\ 

I??2  " * XIj 

If  either  of  these  values  exceeds  the  previous  maximum  for  level  j,  it  is  stored. 
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At  the  end  of  the  simulation,  the  maximum  ?va  lues  of  the  injury  potential 
functions  encountered  are  printed  out  in  graphical  form.  Since  Eq.  (4.11) 
normalizes  IPF  relative  to  1.0  for  each  vertebral  level,  the  values  of  IPF 
indicate  the  likelihood  of  injury  at  each  level.  If  IPF  is  much 'smaller  than 
1.0  at  a level,  injury  is  very  unlikely. whereas  values  of  IPF  in  the  neighborhood 
of  1.0  or  greater  than  1.0  indicate  that  vertebral  failure  is  likely  at  that  level. 

In  order  to  interpret  the  likelihood  of  injury  in  a statistical  sense, 
the  normal  distributions  as  fit  by  Payne  (1972)  to  the  data  of  Geertz,  Crocker 
and  Higgins  and  Nachemson  were  used  to  obtain  standard  deviations.  These  lead 
to  probabilities  of  failure  as  shown  in  Fig.  6.20. 

The  extension  of  this  post-processor  to  other  modes  of  motion  segment 
failure  will  only  involve  the  development  of  additional  Injury  potential 
functions.  Thus,  to  consider  vertebral  dislocations,  it  is  only  necessary  to 
develop  a function  and  appropriate  injury  thresholds  associated  with 
dislocations.  As  part  of  the  function  development,  the  mechanism  of  failure 
and  causative  forces  must  be  determined.  This  procedure  is  currently  under 
study  under  a separate  eohtract  by  the  AMRL. 
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; CHAPTER  V . . 

PRIMATE  HEAD-SPINE  MODELS 

Because  of  the  difficulty  of  obtaining  experimental  results  for  human 
subjects,  data  were  obtained  and  prepared  for  similar  models  of  primates. 

The  goal  of  this  effort  is  to  compare  the  predictions  of  model  simulations 
of  nonhuman  primates  with  experiments  in  order  to  validate  the  basic  modeling 
concepts.  In  addition,  it  is  hoped  that  such  studies  will  shed  further  light 
on  the  procedures  of  scaling  results  from  nonhuman  primates  to  humans. 

In  this  investigation,  preliminary  model  data  was  prepared  for  the 
following  primates: 

1 . chimpanzee ; 

2 . baboon ; 

3.  rhesus  monkey. 

We  will  here  describe  the  sources  of  this  data,  how  it  was  processed,  and  give 
some  results  for  the  models.  In  addition,  we  will  list  some  of  the  major  data 
needed  to  lay  a sounder  groundwork  for  these  models;  it  is  interesting  to 
observe  that  while  data  is  much  easier  to  procure  for  nonhumans,  much  less  is 
reported  in  the  literature. 

V . I Chimpanzee 

Geometric  data  for  the  chimpanzee  spine  was  based  on  Kazarian,  et  al.,  (197$), 
who  provided  the  vertebral  body  heights,  intervertebral  disc  heights  and 
endplate  areas,  and  figures  of  a chimpanzee  spine.  These  figures  were 
digitised  to  provide  the  geometric  data  for  the  model.  Since  only  the  areas 
of  endplates  were  given,  a circular  cross-section  ws?  assumed  and  the  radius 
was  computed.  The  model  structure  is  shown  in  Fig.  5.1. 

The  mass  data  for  the  chimpanzee  was  determined  on  the  basis  of  measurements 
reported  by  Rholes  and  Fineg  (1961),  for  two  groups  of  chimps:  one  group 
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Figure  5.1.  Chimpanzee  spine  model:  (a)  sagittal  plane  view  (b)  frontal  plane  viqw 
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weighing  20-33  lbs.,  the  other  group  weighing  38-54  lbs.  The  goal  of  this 
effort  was  to  obtain.4a.ta  for  s chimpanzee  weighing  100  to  110  lbs,,  so  the 
measurements  in  Rhol&^Fand  Fineg  were  used  only  as  a guide  to  weight  distribution 


Furthermore,  it  had . to^be  assumed  that  weight  distribution  dpes  not  change 

j , . . " - 

significantly  with  change  in  total  body  weight. 

The  procedure  used  to  determine  the  model  masses  is  as  follows.  First 
the  anatomy  of  the  chimps  was  approximated  by  a set  of  geometric  solids  as  shown 
in  Fig.  5.2.  The  volume  of  each  geometric  figure  was  then  computed  using  the 

.. . 3 

measurements  given  in  Rfaoles  and  Fineg,  assuming  a density  of  1 gram/cra  . 

A total  weight  was  computed  for  the  maximum,  minimum,  and  mean  values  of  these 
measurements,  and  a ratio  of  the  weight  of  each  part  of  the  anatomy  to  the  total 
weight.  These  ratios  are  given  in  Table  5.1.  As  can  be  seen,  except  for  a 
few  discrepancies,  the  proportion  between  total  weight  and  the  weight  of  a 
portion  of  the  anatomy  is  somewhat  independent  of  the  total  weight,  so  that 
these  ratios  can  be  used  to  estimate  the  weight  of  the  anatomical  segments. 

Next,  it  was  established  that  for  the  torso,  the  ratio  of  the  width  of  the 
torso  to  the  height  of  the  torso  is  reasonably  independent  of  the  total  weight. 

By  using  a ratio  of  width  of  to  height  of  chest  of  1.18  as  determined  from 
the  measurements  given  in  Rholes  and  Fineg,  the  tor sal  height  of  the  chimp 
for  which  a skeletal  figure  was  provided  by  Kazarian  (1976)  was  estimated. 

Then  using  these  dimensions  and  the  geometry  of  the  torso  as  depicted  in 
Fig.  5.3  (provided  by  Kazarian  (1976)),  the  weight  of  the  torso  segments 
from  Tl  to  L4  were  determined.  Also,  the  total  weight  of  the  chimpanzee  for 
which  the  skeletal  figure  was  provided  was  estimated  and  was  found  to  be  111  lbs. 

To  determine  the  mass  of  each  vertebral  segment  of  the  torso,  if.  was 
necessary  to  know  the  cross-sectional  area  of  the  torso  and  the  thickness  at 
each  level.  The  thickness  of  each  level  was  based  on  the  vertebral  body  heights 
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Figure  5.2.  Idealizations  of  anatomical  components  of  chimpanzee  into  geometric 
solids;  numbered  measur aments  refer  to  Rholes  and  Fineg  (1961). 


Table  5.1  Weight  Ratios 
(each  number  indicates  the  weight  ratio  of  the 
■'  1 ■■  anatomical  component  to  the  total  weight  ) 


4. 

Group  1 ) 

. t 

U 

Group  11 

Med 

Max 

Mia 

Med 

Max 

Min 

Average 

Cheat  - ■ 

0.235 

0.230 

0.264 

0.227 

0.177 

0.237 

0.228 

Lower  torso 

0.166 

0.122 

0.126 

0.145 

0.134 

0.146 

0.143 

Head 

0.137 

0.135 

0.171 

0.112 

0.096 

0.138 

0.132 

Upper  arm 
(each) 

0.033 

0.042 

0.040 

0.039 

0.041 

0.047 

0.040 

Lower  arm 
(each)  • 

0.025 

0.020 

0.007 

0.026 

0.023 

0.036 

0.023 

Hand  (each) 

0.018 

0.042 

0.007 

0.030 

0.067 

0.011 

0.029 

Buttocks 

0.076 

0.060 

0.037 

0.058 

0.082 

0.055 

0.066 

Upper  leg 
(each) 

0.076 

0.071 

0.079 

0.082 

0.071 

0.073 

0.075 

Lower  leg 
(each) 

0.025 

0.022 

0.033 

0.026 

0.022 

0.036 

0.027 

Foot  (each) 

0.014 

0.029 

0.013 

0.026 

0.033 

0.015 

0.022 
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arid  intervertebral  disc  heights  which  were  taken  from  measurements  of  the  skeletal 
figure.  The  area  was  found  frost  Figure  5.2 . From  the  volume  of  the  torso  at 

a vertebral  level,  the  mass  was  computed  by  multiplying  by  the  density  of 

3 • i-:  . 

1 gram/cm  . 

The  equations  for  the  moments  of  inertia,  I,  I , and  I were  derived  on 

x y 2 

the  basis  of  elliptical  cross  sections.  The  results  were  compared  with  the 
masses  of  the  vertebral  levels  of  the  human  spine  used  in  the  Air  Force  model. 

The  masses  of  the  human  spine  were  between  1.2  to  2 times  larger  than  that  for 
the  chimp.  The  human  anatomy  considered  in  the  Air  Force  model  has  a total 
weight  of  169  lbs.,  which  is  approximately  1.6  times  that  of  the  chimpanzee. 

Ihe  masses  and  moments  of  inertias  are  given  in  Table  5.2. 

The  stiffness  data  for  the  intervertebral  discs  was  determined  by  using 
the  procedures  of  Belytschko,  et  al.,  (1974)  and  Schultz,  et  al.,  (1974).  In 
these  procedures,  it  is  assumed  that  the  level-to-level  variations  in  inter- 
vertebral disc  stiffnesses  depend  directly  on  the  geometries  of  the  inter- 
vertebral discs  and  can  be  deduced  from  strength  of  materials  considerations. 
Furthermore,  since  no  disc  stiffnesses  were  available  for  chimpanzees , the 
material  properties  were  assumed  to  be  identical  to  human  discs.  The  axial, 
bending  and  torsional  stiffnesses  at  a level  J are  given  by 
(axial  stiffness)^  * 

B 

(bending  stiffness)^  » R^  Kj  (5.1) 

(torsional  stiffness)^  * 


where 


A 


K 


J 

5 

J 


IDH 


0 
IDH 


5 / 4 
l r 
\ o 


2 K 


B 


(5.2) 
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Table  5.2b.  Inertial  data  for  vertebral  levels  L4-TU  for  chimpanzee 
with  viscera  model 


Level 

”<2)  | 

I~(3) 

t.<3) 

y 

- (3V 

z 

Visceral 

nasa  ds*-a 

L4-T11 

L4 

6.690 

1.150 

3.215 

4.281 

L3 

8.870 

2.042 

4.483 

6.394 

L2 

10.030 

2.581 

5.431 

7.874 

T.l 

11.570 

3.214 

6.785 

9.841 

112 

1.0.820 

2.782 

6.340 

8.975 

Til 

8.100 

1.855 

4.722 

6.515 

Vertebral 

nass  data 

L4-T11 

L4 

2.310 

0.128 

0.239 

0.338 

L3 

2.230. 

0.091 

0.267 

0.325 

L2 

2.060 

0.067 

0.292 

0.331 

Ll 

1.960 

0.060 

0.273 

0.306 

T12 

2.710 

O.OSi 

0.717 

0.762 

Til 

| 2.030 

0.055 

0.531 

0.57). 

■ -a 

Inertial  data  tor  torso  segments  T10-L3  is  given  in  Table  5.2a 

2 

Translational  mass  - grass  x 10 

- - - 2 

Rotational  masses  about  coordinates  x,  y,  z gta  cs  x 10 


,/s*; 


8fe 


m.. 


>.V#; 


where  r.  and  r are  the  inside  and  outside  radii  of  the  annulus  fibrosis, 
i o * 

and  IDH  is  the  intervertebral  disc. height,  Eqs.  (5,2)  fellow  from  elementary 

strength  of  materials  treatment  of  the  stiffness  of  a disc,  it  was  assumed 

that  r,  * 0.75  r . 
i o 

The  geometric  factors  of  the  chimp  spine  were  based  on  data  provided  by 
Kazarian,  We  assumed  that  the  R.^,  Rp,  and  R^  for  chimpanzees  and  humans  were 
the  same,  i.e,  that  material  properties  of  the  two  species  are  similar.  These 
constants  were  then  determined  in  the  following  manner.  In  Table  1 of  Belytschko, 
et  al,,  (1974),  geometric  data  at  intervertebral  levels  L2/L3,  T12/L1,  TlO/Tll, 
and  T8/T9  were  given.  From  this  data  A,  and  for  those  levels  were 
calculated. 

Table  2 of  Schultz,  et  al.  gives  stiffnesses  at  T,2/T3,  T12/LI,  TlO/Tll, 
and  T8/T9,  Ratios  were  set  up  for  corresponding  quantity ; c the  corresponding 
levels  and  the  constants  R^,  R^  and  R^  were  determined  using  the  average  of  the 
four  levels  (corresponding  values  at  all  four  levels  deviated  little  from  the 

respective  averages).  The  values  are 

9 2 

= 0.144  x 10  dynes/cm 

Rp  ® 0.746  x 10^  dynes /cm^  (5.3) 

8 2 

“ 1.27  x 10  dynes /cm 

T B 

The  geometric  data  (A,  K,  K ) for  the  chimp  appears  in  Table  5.3  while  the 
stiffnesses  for  the  chimp  spine  model  appear  in  Table  5.4, 

It  should  be  noted  that  there  were  discrepancies  between  the  tabulated 
geometry  (used  in  obtaining  the  stiffnesses)  and  the  geomtrey  of  the  chimpanzee 
skeletal  figure.  The  reasons  tor  this  are: 

1)  IDH  for  the  giv' n data  is  larger  than  the  IDH  measured  from  the  chimp's 
skeletal  figure.  There  is  about  a 50%  discrepancy. 


% 
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Table  5.3 


Geometries  of  Chimpanzee  Intervertebral  Piece 


avel 

IDH (cm) 

2 

Area (cm  ),A 

r (cm) 

0 

ri(cm)  ! 

C3 

0.42 

2.94 

0.97 

0.73 

C4 

0.42 

2.87 

0.96 

0.72 

C5 

0.42 

2.89 

0.96 

0.72 

C6 

0.43 

2.66 

0.92 

0.69 

C7 

0.40 

2.74 

0.93 

0.70 

Tl 

0.71 

3.10 

0.99 

0.74 

T2 

0.73 

3.28 

1.02 

0.77 

T3 

074 

3.37 

1.04 

0.78 

T4 

0.74 

3.64 

1.08 

0.81 

T5 

0.74 

3.43 

1.04 

0.78 

T5 

0.70 

4.86 

1.24 

0.93 

T7 

0.94 

4.31 

1.17 

0.88 

T3 

0.94 

4.22 

1.16 

0.87 

T9 

0.95 

5.05 

1.27 

0.95 

T10 

1.00 

5,73 

1.35 

1.01 

Til 

1.00 

5.97 

1.38 

1.04 

T 12 

1.20 

6.44 

1.43 

1,07 

Ll 

i.30 

7.29 

1.52 

1.14 

L2 

1.40 

8.05 

1.60 

1.20 

U 

1.20 

8.10 

1.61 

1.21 

L4 

1.40 

9.43 

1.73 

1.30 

hare  IDH, 

Area  are 

given 

ro  = J 

Area 

rt 

B 1 

K " 2 IDH 

, 4 4, 

<Co  ' rl> 

1,4  4. 

IDH  (ro  - ri) 


3 

4 ro 
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Table  3,4 

Chimpanzee  Stiffness  Data 


Intervertebral 


Discs 


Visceral' 

Elements 


Refers  to  element  below  designated  vertebral  level 


Visceral  elements  are  located  between  the  pelvis  and  TlO.  Stiffnesses  of 
all  viscera/spine  interconnecting  elements  « 1.0  x 10^  dyne/cta 


Though  based  on  data,  appears  spurious 
1 kp  w 106  dynes 
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2)  Area  for  the  given  data  is  smaller  than  the  area  measured  from  the  chimp 
skeletal  figure.  There  is  about  a 25%  discrepancy. 

The  chimp  stiffnesses  in  the  upper  spine  are  only  0-4  of  the  human  stiffnesses. 
However,  the  chimp  data  lacks  the  relative  maximum  stiffness  found  in  humans 
around  Til,  but  instead  increases  continuously  to  L4,  where  the  stiffnesses 
are  comparable  to  human  intervertebral  disc  stiffnesses. 

The  visceral  bulk  modulus  was  obtained  by  using  Eqs.  (4.5)  and  (4.6)  of 
Belytschko,  et  al.,  (1976),  using  a radius  of  12  cm  as  taken  from  Table  5.2 
and  a membrane  thickness  of  0.5  cm.  Tee  resulting  stiffnesses  are  given  in 
Table  5.4. 

A modal  analysis,  an  impedance  analysis  and  an  ejection  analysis  of  the 
model  of  a chimpanzee  spine  were  performed  with  the  ligamentous  spine  both 
with  and  without  viscera.  Table  5.5  lists  the  first  six  axial  natural  frequencies 
and  the  first  seven  sagittal  plane  natural  frequencies  for  the  chimpanzee  isolated 
ligamentous  spine  and  compares  them  with  the  AMRL-TR-76-10  human  spine  model. 
Although  a casual  comparison  indicates  that  the  frequencies  of  the  chimpanzee 
spine  in  the  axial  mode  ere  somewhat  smaller  than  those  for  the  model  of  the 
human  spine,  a comparison  of  the  modes  (not  given)  shows  that  there  is  a 
significant  difference  in  the  nature  of  the  lower  modes  between  the  chimpanzee  and 
the  human  models.  In  the  chimpanzee,  the  two  lowest  axial  modes  are  almost 
entirely  associated  with  head  movement  and  involve  very  little  deformation  of 
the  lover  thoracic  and  lumbar  regions,  whereas  in  the  human  spine  model  the 
lowest  axial  mode  involves  considerable  lower  back  deformation.  Thus  the 
chimpanzee  axial  mode  corresponding  tc-  the  first  human  axial  mode  of  17  Hz 
is  54  Hz.  Thus,  as  expected,  the  chimpanzee  racdel  modes  are  significantly 
higher.  Similarly,  the  first  flexural  modes  of  the  chimpanzee  involve  primarily 
head -neck  motion. 
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TABLE  5,5 


Comparison  of  Natural  Frequencies  of 


Human  Spine  Model  (HSM)  and  Chimpanzee  Spine  Model  (CSM) 


Axial  Frequencies  (Hz') 

HSM 

CSM 

17.09 

13.97 

32.31 

31.81 

51.29 

53.96 

77.22 

77.72 

100.74 

105.22 

124.98 

131.91 

Sagittal  Plane  Frequencies(Hz) 

HSM 

CSM 

Ratio 

/CSM\ 

\HSM/ 

1.28 

0.76 

0.59 

3.14 

2.04 

0.S5 

5.99 

3.70 

0.62 

9.94 

5.32 

0.54 

13.31 

7.24 

0.54 

16.71 

8.55 

0.51 

18.45 

10.30 

0.56 
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This  difference  between  the  chimpanzee  and  human  can  be  directly  related 
to  the  mass  distribution.  In  the  chimpanzee,  the  ratio  of  head  mass  to  total 
mass  is  quite  large  compared  to  the  human,  whereas  the  opposite  holds  for  the 
pelvic  mass.  It  should  be  noted  that  these  differences  in  mass  in  our  model 
are  based  on  rather  crude  data  obtained  from  gross  measurements.  However, 
this  difference,  while  not  exactly  represented  in  our  model,  certainly  exists, 
and  it  could  shed  some  light  on  the  difficulties  in  scaling  studies;  because 
of  the  differences  in  mass  distributions,  the  modal  character  of  the  human 
and  chimpanzee  are  evidently  significantly  different,  thus  precluding  any 
simple  scaling  laws. 

An  ejection  simulation  was  made  with  the  chimpanzee  spine  - viscera  model. 
The  prescribed  acceleration  consisted  of  a peak  magnitude  of  lOg  which  was 
reached  in  14  msec,  with  an  onset  rate  of  714g  per  second.  The  lOg  magnitude 
was  maintained  at  constant  value  for  the  remainder  of  the  simulation.  The 
total  duration  of  the  simulation  was  40  msec.  Figure  5.4  shows  the  axial  force 
time  histories  at  T9-T10  and  T3-T4. 

Impedance  curves  were  obtained  for  the  two  column  chimpanzee  model.  These 
are  shown  in  Fig.  5.5. 

We  would  again  like  to  stress  the  large  number  of  extrapolations  that  were 
necessary  in  obtaining  the  100  lb.  chimpanzee  model.  Any  additional  data 
would  be  of  great  value.  Some  suggested  data  are; 

a)  axial  stiffness  at  any  level  for  a chimpanzee  cadaver  motion  segment 
(it  would  be  useful  to  have  the  size  or  weight  of  the  chimp  corresponding  to 
this  stiffness); 

b)  any  other  disc  stiffness,  e.g.  bending; 

c)  a series  of  photographs  or  x-rays  of  a 100  lb  chimp  to  verify  our 
size  scaling  procedures. 
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Time  (sec) 


Figure  5. 


. Axial  force  response  in  the  T3-T4  and  S^-TlG  discs  from  ejection 
of  chimpanzee  model. 
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MODULUS  OF  IMPEDANCE 
■ DYNE-SEC/CM  x 10s 


20  40  60  80  100  rad/sec 

3.2  6.4  9.5  ! 2.7  15.9  Hz 


FREQUENCY 

Figure  5.5  Impedance  of  chimpanzee  modal  with  viscera:  (1)  pelvic  mass  impedance; 
(2)  impedance  of  chimpanzee  without  pelvic  mass;  (3)  impedance  of 
chimpanzee  with  pelvic  mass;  (b)  impedance  of  chimpanzee  with  buttocks. 

7 

Buttock  stiffness  is  6.55  x 10  dyne/cm. 


V.2  Baboon 

The  geometrical  data  of  the  baboon  was  obtained  in  the  same  manner  as 
that  of  the  chimpanzee,  based  on  data  of  Kazarian  (1976).  The  resulting 
model  geometry  is  plotted  in  Fig.  5.6. 

Mass  data  is  baaed  on  Reynolds  (1974),  who  measured  the  inertial  properties 
of  anatomical  segments  of  four  frozen  cadaver  baboons.  The  baboons  ranged  in 
mass  from  10.8  kg  to  12,66  kg.  It  is  of  interest  to  observe  that  the  densities 
were  also  measured  by  Reynolds  (p.  73);  they  ranged  from  a low  of  0.?S  gm/cm 

3 

in  the  torso  to  a high  of  1.12  gm/cm  in  the  feet  and  hands,  with  an  average 

3 

of  1.016  gm/cm  , Thus  the  density  assumptions  made  in  determining  the  chimpanzee 
inertias  are  quite  reasonable:  the  major  sources  of  error  would  probably  be 
the  idealizations  of  anatomical  segment  volumes  and  the  extrapolations.  The 
torso  percentage  of  total  volume  in  the  baboon  is  higher  than  our  estimated 
percentage  for  the  chimpanzee  (537.  vs  447.)  and  the  head  is  lower  (8.37,  vs  13.27). 

The  head  mass  and  moments  of  inertia  were  taken  directly  from  the 
measurements  of  Reynolds.  Since  Reynolds  did  not  segment  the  torso,  the  inertial 
properties  of  the  torso  segments  were  obtained  by  using  the  data  of  Kazarian 
for  vertebral  body  and  intervertebral  disc  heights,  a radius  of  7.35  cm  and 
the  density  given  by  Reynolds.  The  inertia  data  is  given  in  Table  5.6. 

Stiffness  data  was  constructed  using  the  measurements  of  Kazarian  (1976), 
Eqs.  (5.1)  and  (5.2),  and  the  R_^  factors  from  human  motion  segments.  This 
data  is  given  in  Table  5.7.  No  separate  visceral  column  is  included. 

The  impedance  of  the  baboon  model  is  shown  in  Fig.  5.7.  It  should  be 
noted  that  here  reliable  mass  data  is  available  for  the  head,  yet  the  model 
exhibits  a substantially  different  impedance  curve  from  the  human  model, 
indicating  again  a potential  source  of  difficulty  in  scaling. 
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TABLE  3.6 
INERTIA  DATA 

2 3 

Mooents  of  Inertia  (go-co  x 10  ) 


•90 

20.240 

’08 

2.708 

'80 

(33 

3.333 

:62 

3.262 

12.37$ 


Pelvis 


16.430 


37.2GO 


11.726 

59.173 


Neck 

0.25 

Tl 

1.90 

12 

1.60 

T3 

1.70 

T4 

1.80 

T5 

2.10 

T6 

2.00 

T7 

2.20 

T8 

2.40 

T9 

2.50 

TlO 

2.80 

Til 

2.90 

Tl2 

3.10 

Ll 

3.00 

u 

4.40 

L3 

4.70 

14 

5.20 

1.5 

5.70 

16 

5.60 

17 

5.90 

TABLE  5.7 


STIFFNESSES  OF  BABOON  DISCS 


ffnees 


Bending 

Q 

( dyne  -cgi  /rad  x 10  i 


Tora tonal 


1.76 

1.56 

1.60 

1.45 

1.07 

0.98 

0.42 

0.51 

0,4i 

0.43 

0.34 

0.34 

0.30 

0.24 

0.28 

0.22 

0.18 

0.18 

0.22 


1.50 

1.33 

1.37 

1.23 

0.91 

0.84 

0.36 

0.43 

0.38 

0.37 

0.29 

0.29 

0.25 

0.20 

0.24 

0.19 

0.15 

0.15 

0.19 


V.3  Rhesufl  Monkey 


For  thi  rhesus  monkey,  only  the  measurements  and  skeletal  figures  were 
available  for  data  extraction.  The  geometry  of  the  model  is  snown  in  Fig.  5.8, 

For  the  rhesus  monkey  a total  weight  of  21  lbs  (9.55  kg  mess)  was  specified,  so 
the  mass  distribution  data  was  scaled  from  the  baboon.  Stiffness  data  was  obtained 
using  the  ame  procedure  -a  for  the  chimpanzee  and  baboon.  Inertial  and  stiffness 
data  is  given  in  Tables  5.8  and  5.9.  The  impedance  of  the  rhesus  monkey  model  is 
a Iso  shov’i  in  i . . 5.9, 

In  addition  to  . fcter  stiffness  deta,  which  is  needed  for  all  of  the  primate 
models,  specific  inertia  data  would  be  useful  for  the  rhesus  monkey. 
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Figure  *.8.  KheSua  acokay  dfcine  ssodei;  (a)  sagittal  plau«  view; 
(b>  frccwl  plaa^  view 
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Mass  m 


Table  5.8,  Rhesus  monkey  inertia  data 

2 3 

Moments  of  Inertia  (gm-cm  x 10  ) 


Table  5,9,  Rhesus  monkey  stiffness  data 
Axial  Stiffness  Bending  Torsional 
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Figure  5.9, 


Icpcdance  of  rhesus  turnkey: 
pelvis;  (3)  with  pelvis 


(1)  pelvic  mass  impedance:  (2)  without 
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CHAPTER  VI 


SIMUIATICNS  OF  PRE -EJECTION  PILOT  ALIGNMENT 
VI. 1 Motivation  and  Method 

Mohr,  et  al.,  (1969)  and  Li,  et  al.,  (1970),  have  discussed  the  importance  of 
the  pre-ejection  posture  as  a factor  in  modifying  a crewman’s  tolerance  to  acceler- 
ations experienced  during  ejection  from  disabled  high  performance  aircraft. 
According  to  Bosee  and  Payne  (1961),  the  optimal  spinal  alignment  is  the  normal 
seated  position,  with  the  vertebral  bodies  “located  squarely  over  each  other".  Any 
deviation  from  this  optimal  spinal  configuration,  such  as  that  resulting  when  a 
crewman  is  slumped  forward  or  sidewards  in  the  seat,  would  increase  the  possibility 
of  injury  during  ejection. 

In  this  Chapter  we  will  .onsider  the  pilot  alignment  phase  of  the  ejection 
sequence  and  will  present  a method  for  analyzing  the  spinal  response  to  loadings 
encountered  during  alignment. 

VI. 2 Computer  Modeling  of  Belt  and  Shoulder 

The  pre-ejection  pilot  alignment  is  accomplished  via  the  powered  inertia  lock 
reel  which,  when  required,  activates  the  restraining  belt,  forcing  the  crew  member 
into  the  eject  position  (Fig.  6.1).  In  the  pre-ejection  alignment  models  the  belt 
forces  are  transmitted  to  the  spine  through  rigid  bodies  representing  the  shoulders 
and  through  the  ribs.  Based  on  the  work  of  Eyeleshyoer  and  Schoeeaker  (1970) , it 
was  determined  that  the  sagittal  plane  cross-sections  of  the  shoulders  could  be 
approximated  as  circular  segments. 

The  belt  forces  are  transmitted  to  the  shoulder  rigid  bodies  and  the  ribs 
through  points  (which  we  shall  call  contact  points)  which  approximate  the  shape  of 
the  upper  torso  in  planes  parallel  to  the  sagittal  (ys)  plane,  8.13  ca  on  either 
aide  of  the  sagittal  plane.  Fig.  6.2  shows  a typical  contact  point  coof iguretion 
in  s plane  parallel  to  the  sagittal  plane.  Contact  points  1 to  j + 2 are  secondary 
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Figure  6.Z.  Typical  contact  print  configuration 
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nodes  associated  with  the  shoulder  rigid  body  primary  node,  called  S&B.  Contact 
points  j + 3,  to  n are  secondary  nodes  associated  vith  ribs,  Ri  to  R6,  respectively. 
Each  rib  has  only  one  contact  point.  The  inertia  lock  reel  (ILR)  and  tie  down 
point  (TOP)  locations  shown  in  Fig.  6.2  are  not  drawn  to  scale. 

At  each  time  step,  the  location  of  a contact  point,  J,  is  given  by 

*kJ(V  “ xk.T<’ti-l)  + ukJ(ti-l)  (6*1) 

where 

*u,  x2j,  x^  are  the  global  coordinates  j (x^  and  x,  and  y and  x3  and 
z are  used  interchangeably) 

u^jCt^)  = displacements  of  contact  point  J computed  at  time  step 

Next  a check  is  made  for  contact  between  the  belt  and  the  contact  points.  Contact 
between  a posterior  contact  point  (toward  the  seatback)  and  the  belt  is  determined 
by  checking  whether  contact  at  that  point  would  introduce  any  concave  kinks  into 
the  belt.  At  all  times,  the  position  of  the  inertia  lock  reel  (ILR)  is  given,  so 
this  is  determined  by  checking  the  sign  of  the  global  x component  of  the  cross 
product  of  the  vector  from  the  inertia  lock  reel  to  a possible  contact  point  and 
:'r:.  vector  from  the  inertia  lock  revl  to  the  nert  contact  point,  as  illustrated  for 
to.itact  point  l in  Fig.  6,2,  where  the  condition  would  be 

:»  0,  belt  does  not  contact  point  1 

(6.2) 

> 0,  belt  contacts  point  1. 

Contact  between  an  anterior  (chest  side)  contact  point  snd  the  bait  is  similarly 
determined  so  that  no  concave  kinks  exist  in  the  on  It  between  the  chest  and  the  tie 
down  point.  For  axampie,  for  contact  point  t»  in  Fig'<  6.2, 

2 0,  belt  does  not  contact  point  n 

(6.3) 

< D,  belt  ccr-accs  point  a. 

No  force  is  applied  to  any  point  which  is  i*ot  in  contact  with  the  belt,  so  loss  of 


l«(vl  X V,)  | 


i.(v3  x v4)  { 
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contact  between  the  belts  end  parte  of  the  f’  juldor  end  torso  la  possible. 

the  force  (in  the  plane  of  the  contact  points)  on  a typical  contact  point,  for 
example  point  J in  Fig#  6.2,  18  determined  by 

?j  * 1 * ;J+i.j>  <«•*> 

where  T = belt  tension,  and  e^  and  ^ are  unit  vectors  along  the  belt 
shown  in  Fig.  6.2. 

The  global  components  of  Fj  are  then  given  by 


aa 


where 


F - t! 

xj+i,j\ 

xj 

aj-i 

lj  } 

F . - t| 

<yj-i.j 

yj  ' 

^ lj-i 

F » T| 

2j+i.j\ 

jy  1 

v Lj-i 

lj  ' 

v 5 

J-1,J 

* xj-i  " xj*  etc- 

<6. 5) 


Lj  " distance  between  contact  points  J+l  and  J. 


Once  the  global  components  of  the  contact  forces  have  been  computed  at  the  contact 
points  they  err  formed  into  forces  and  moments  acting  at  the  primary  nodes 

(mass  centers)  of  the  shoulder  rigid  bodies  and  the  ribs. 


The  resultant  force  at  « shoulder  rigid  body  primary  node  is  given  by 
m 


RF 


j 


1-1 


n 


(6.6) 


where  ts  * number  of  contact  points  defining  the  Shoulder  rigid  body.  The  rcsulc&ot 


moment  at  « shoulder  rigid  body  primary  node  is  given  by 
e 


{«)■-[  e COj]  t'.f  frjl] 


(6.6) 


woere 


[Oil 


1 0 

51 

•Ft 


•SI  yi! 


xi  0 


(6.8) 
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and  arc  expressed  in  the  body  coordinate  system  of  the  coordinated  of  primary 

node  P,  and  [Xl  is  the  body  to  global  coordinate  system  transformation  matrix  for 
the  shoulder  rigid  body  primary  node  P.  Since  each  rib  has  only  one  contact  point, 
the  forces  and  moments  at  the  rib  primary  nodes  are  given  by  Eqs.  (6.6)  and  (6.7) 
without  the  summations. 

The  belt  forces  on  the  contact  points,  once  they  are  transformed  into  forces 
and  moments  at  the  primary  nodes  of  the  shoulder  rigid  bodies,  are  treated  like  the 
other  external  forces  in  the  equations  of  motion.  No  predictor  programming  is 
included,  so  the  belt  model  is  only  suited  for  explicit  time  integration. 

The  morphology  of  the  interconnections  of  shoulders  and  ppine  is  rather  compli- 
cated. Furthermore,  data  on  the  various  components  of  the  inter connective  tissues 
are  not  available.  Since  forces  in  these  in ter connectors  themselves  were  not  of 
Interest,  they  were  simplified  as  much  as  possible. 

The  shoulder  rigid  bodies  were  connected  by  beam  elements  to  Tl,  T2  and  T3,(see 
Fig.  6.3),  and  the  stiffnesses  of  the  beam  elements  were  chosen  on  the  basis  of  a 
parameter  study.  The  stiffnesses  were  varied  so  that  when  the  shoulders  are  pulled 
back,  the  spinal  displacement  would  lag  behind  the  shoulder  displacements  by  a small 
amount,  but  so  that  no  snapback  or  overshoot  in  the  spinal  response  was  observed. 

The  stiffnesses  of  these  elements  and  locations  of  the  contact  points  are  given 
in  Table  6.1. 

VI. 3 Pre-Ejection  Pilot  Alignment  Models 

In  the  following,  we  will  describe  how  the  spine  models  were  prepared  for  the 
simulation  of  the  pre-ejection  alignment  runs.  Ue  will  designate  these  models  as 
PAMX-Y,  where  PAW  denotes  Pre-ejection  Alignment  Model  and  X indicates  the  initial 
forward  displacement  of  Tl,  and  Y (if  present)  indicates  an  initial  lateral  displace* 
menu  The  models,  PAH2,  ?AH4  and  FAH4-4,  were  constructed  from  the  Isolated,  liga- 
mentous spine  without  c rib  cage.  In  2AM6,  the  entire  model  was  included.  PAM2 
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Table 

6.1,  Coordinates  of  contact  points  and  stiffnesses  of 
connecting  the  shoulder  rigid  bodies  to  vertebra 

beam  elements 
Tl,  T2  and  T3 

LSRB^  primary  node  «lobal  coordinates  (cm): 

x = 8 

.130,  y ® 

-13.943, 

z « 42.900 

Local 

(21 

coordinates  of  contact  points  (CP) 

C? 

A 

X 

y 

A 

Z 

1 

0.0 

6.447 

2.725 

2 

0.0 

4.221 

5.533 

3 

0.0 

0.864 

6.946 

^r- 

4 

0.0 

-2.726 

6.448 

5 

0.0 

-5.584 

4.222 

6 

0.0 

-6.947 

0.865 

7 

0.G 

-6.448 

-2.725 

Global  coordinates  of  left  rib  p;  !• 

•<vary  nodes  and  local 

coordinates  1 

of  associated  contact  point  (CP) 

Rib 

X 

y 

z 

CP  x 

A 

y 

A 

Z 

t 

•4 

2.556 

-16.328 

41.569 

8 5.574 

-2.871 

-5.069 

2 

3.380 

-16.355 

37,915 

9 4.750 

-2.245 

-5.515 

3 

4.115 

-15.548 

34.556 

10  4.015 

-2.938 

-6.374 

4 

5.121 

-14.320 

31.460 

11  3.009 

-3,827 

-7.033 

5 

5.548 

-14.089 

27.768 

12  2.782 

-3.535 

-6.305 

6 

5.574 

-12.958 

24.47S 

13  2.556 

-3.960 

-6.255 

Stiffnesses  of  beats  eles«nts  between  shoulder  rigid  bodies  and  TI-T3 


Axial : 
Torsion: 
Bendi ng : 


5 x 10  dyne-co 
7 

l x 10  dyue-ca 

kl  * M2 

k.j  e l x 107  dyne  -ers 
k^  «*  1 x 1010  dyae-ce/racT 


l^LSRB  • left  shoulder  rigid  body 
(2> 

Origin  located  at  LSRB  pr Leary  node 
^Origin  located  at  primary  node  o£  rib 

The  right  hand  side  configuration  is  s oirror  iaage  of  the  left  side. 
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(Fig.  6.4)  was  obtained  by  taking  the  isolated  thoraco-lumbar  spine  model (Fig.  2.1X 
and  placing  Tl  5.08  cm  (2  in)  forward  from  its  undeft-rmed  position,  which 
is  approximately  9.3  cm  (3.7  in)  forward  from  the  seatback.  The  remaining  vertebral 
bodies  were  then  positioned  so  as  to  obtain  a smooth  curve  for  the  spine,  which  is 
shown  in  Fig.  6.4. 

?AM4  (Fig.  6.5)  was  obtained  from  the  normal  configuration  of  the  isolated 
thoraco- lumbar  spine  by  a static  analysis  in  which  Tl  was  displaced  forward  from 
its  undeformed  position  by  10.16  cm.  (4  in).  The  static  analysis  then  determined 
the  remaining  vertebral  positions  by  equilibrium.  The  initial  location  of  Tl  is 
approximately  14.4  cm  (5.7  in)  from  the  seatback. 

In  PAM4-4  (Figs.  6.6  and  6.7)  vertebra  Tl  was  displaced  4 inches  forward  and 
4 inches  to  the  left.  The  initial  configuration  was  obtained  in  the  same  manner 
as  PAM4. 

PAMo  (Fig.  6.8)  was  obtained  by  taking  the  complex  spine  model,  CSM  (Fig.  2.2) 
and  rotating  it  forward  about  L5  until  Tl  was  15.24  cm.  (6  in)  forward  from  its 
undeformed  position.  This  is  approximately  19,4  ere  (7.7  cm)  from  the  seatback. 

VI, 4.  Pre-Election  Pilot  Alignment  Studies 
VI. 4.1  Study  1. 

In  the  first  study,  the  shoulder  rigid  bodies  were  modeled  with  6 contact  points 
per  shoulder,  as  shown  in  Fig.  6.9  and  6.10.  The  shoulders  are  aligned  symmetrically 
with  respect  to  the  seatback.  The  reel  is  located  at  approximately  the  standard 
height,  3.4  cm  (1.3  in)  above  the  acromion.  From  Kycieshymer  and  Schoemaker (1970) , 
the  acromion  appears  to  be  approximately  3.2  cm  (1.3  in)  above  the  centroid  of  Tl, 
while  from  Pansky  and  House  (1964),  the  acromion  appears  to  coincide  with  the  bottom 
of  Tl . The  acromion  height  in  the  spine  models  was  taken  to  be  the  average  of  these 
two  values.  The  tie  down  points  are  8.13  cm  (3.2  in)  on  either  side  of  the  pelvic 
mass  center.  The  force  in  each  belt  was  aero  initially  and  increased  linearly  to 
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Response  of  PAM2 


2.22  x 10®  dynes  (500  lbs)  in  150  msec,  for  a combined  belt  force  of  4.45  x 10® 
dynes  (1000  lbs)  at  150  msec. 

The  models  PAM2  and  PAM4  were  used  in  this  study.  Figures  6.9  and  6.10  show 
the  0,  50,  100  and  150  msec  pilot  configurations  for  PAM2  and  PAM4  respectively. 
Table  6.2  compares  peak  axial  force  and  sagittal  plane  bending  moment  magnitudes 
for  the  neck  and  intervertebral  discs  T3/T4  and  L3/L4. 

VI. 4, 2 Study  2. 

Study  2 considers  a nonsymmetric  initial  configuration.  PAM4-4  was  subjected 
to  the  same  pre-ejectiou  belt  force  as  that  in  Study  1.  Figures  6.11  and  6.12 
depict  the  0,  50,  100  and  150  msec  pilot  configurations,  sagittal  and  frontal  plane 
views  respective!/.  Two  features  of  the  model  had  adverse  effects  on  the  results: 

(1)  the  orientations  of  the  shoulder  rigid  bodies  were  rotated  too  much  initially; 

(2)  the  seatback  does  not  offer  any  resistance  to  lateral  motion  of  the  spine 
because  Coulomb  frictional  resistance  is  not  modeled. 

VI. 4. 3 Study  3. 

Here  the  effects  of  changing  the  location  Of  the  reel  were  studied.  Two  runs 
were  made,  one  with  the  reel  2.2  cm  (0.9  in)  below  the  standard  height,  and  one 
with  the  reel  7.2  cm  (2.8  in)  below  the  standard  height.  The  model  used  was  PAM2 
with  the  mass  eenter  of  the  head  shifted  slightly  closer  to  the  seatback  than  in 
Study  1,  The  tie  down  point  location,  belt  force  time  history,  and  shoulder  rigid 
bodies  orientations  were  the  same  as  in  Study  1.  Figures  6.13  and  6.14  show  the  0, 
SO,  100  and  150  msec  pilot  configurations  for  the  reel  at  50, § «a  and  45.8  cm 

respectively, 

VI. 4. 4 Study  4. 

Here  the  effects  of  removal  of  a portion  of  the  seatback  was  considered.  The 
model  used  was  PAM2,  The  reel  height  was  0.9  in  below  the  standard  height  and  the 
tie  dovr.  point  location,  belt  force  cia 2 history  and  shoulder  rigid  bodies  orien- 


Flgur«  6.11.  Study  2:  Sagittal  plane  response  of  PAM4-4„ 


Piguru  6,12,  Study  2:  Front*!  pienc  re«pon*e  of  PAK4-4. 


Plgwre  6 .13.  Study  3:  Responae  of  PAM2,  reel  height  2.2  cm  below  standard 
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tation  were  the  same  as  in  Study  1.  The  aeatback  was  "cut-off"  betweeen  T5  and  T6 
so  that  the  head  and  vertebral  bodies  T1-T5  were  free,  Figure  6,15  shows  the  0,  50, 

100  and  150  m3ec  pilot  configuration*.  

VI. 4. 5 Study  5.  , 

In  study  5,  PAM4  was  subjected  to  the  belt  force  time  history  shown  in  Fig.  6.16. 
The  continuous  curve  is  that  of  AMRL  retraction  test  run  No.  93  for  which  the  retrac- 
tion distance  was  6 in.  The  reel  is  locited  at  the  standard  height  ar i the  tie  down 
point  location  is  the  same  a?  in  the  previous  studies.  The  shoulder  rigid  bodies, 
which  are  aligned  symmetrically  with  respect  to  the  sestback,  have  been  refined  to 
more  closely  approximate  the  shape  of  the  shoulders  and  are  modeled  by  7 contact 
points  per  shouluer  rigid  body.  Although  the  shape  of  the  shoulder  rigid  bodies  is 
approximated  as  a circular  .1egme.1t,  this  shoulder  rigid  body  is  more  refined  than 
those  u.*ed  in  the  previous  su.aies  and  the  radius  was  chosen  to  agree  with  the 
actual  physical  dimensions.  Figure  6,17  depicts  the  0,  50,  100,  150  and  200  msec 
pilot  configurations.  Table  C.3  lists  the  peak  axial  compressive  forces  ana  sagittal 
plane  bending  moments  and  the  approximate  times  at  which  they  occur. 

VI. 4.0  Study  6, 

In  the  final  study,  the  rib  cage  was  included,  and  belt  interaction  with  both 
the  shoulder  rigid  bodies  and  the  ribs  was  considered  to  investigate  the  effects  of 
belt  force  transmission  through  the  rib  cage.  The  reel  and  tie  down  point  locations 
were  the  same  as  for  study  5.  The  belt  force  time  history  is  again  Ch&t  shown  in 
Fig.  6.16.  Figure  6.18  is  a photograph  of  an  initial  pilot  configuration  used  for  the 
AMRL  retraction  test  run  No.  93  for  which  the  retraction  distance  was  6 in.  The 
initial  pilot  configuration  used  in  this  study,  PAM6  (Tl  6 in  forward  from  its 
normal  position,  7.7  in  forward  from  the  seatback),  eppear&  to  be  similar  tc  that 
shown  in  Fig.  6.18.  but  the  photograph  only  gives  a rough  idea  of  the  spine's 
initial  configuration.  Figure  6.19  shows  the  0,  50,  100,  150  and  200  msec  pilot 


Table  6.3  Maximum  internal  forces  for 

pre-ejection  alignment  study  5. 


Level 

Maximum 

force 

F 

Time  of 
max.  force 
t(F> 

Moment 

at 

. t(F) 

Maximum 

moment 

M 

Time  of 
max.  moment 
t(M) 

Force 

at 

t(K> 

L5-L4 

16.13 

128 

2 

7.90 

171 

8 

L4 

15.76 

126 

5 

9.76 

142 

4 

L3 

16.18 

130 

6 

10.38 

138 

10 

L2 

16.77 

130 

6 

8.23 

112 

4 

I.l 

16.66 

130 

5 

8.48 

110 

1 

T12 

14.13 

128 

3 

7.96 

110 

i 

Til 

13.13 

143 

6 

7.77 

169 

2 

TlO 

11.56 

143 

1 

9.34 

108 

1 

T9 

10.83 

140 

4 

8.86 

185 

8 

T8 

9.74 

140 

10 

9.55 

138 

9 

T7 

10.92 

138 

7 

9.48 

132 

6 

T6 

11.06 

124 

8 
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Figure  6.18.  Initial  pilot  configuration  for  AMRL  retraction  test  run  No.  93. 
Retraction  distance  is  6". 


Figure  6.19.  Study  6:  Response  of'PAjJfr 


Before  considering  any of ‘the ^il^'^lignaent  studies  individually, 
some  coranents  on  certain  aspects  of  modeling  should  be  made.  In  all  of  these 
studies,  the  neck  is  modeled  as  a single  beam  element  connecting  the  head  directly 
to  Tl.  The  flexibility  of  this  neck  element  results  in  exaggerated  motions  of  the 
head,  leading  to  excessive  rotation  of  Tl.  This  is. particularly  evident  for  the 
cases  with  Tl  more  than  2 in  forward  from  Jits  normal  position  (e.g.  studies  5 and  6). 
The  muscles  of  the  neck  would  distribute  the  head  inertial  effects  over  a larger 
part  of  the  torso,  thus  reducing  the  moments  at  Tl.  Furthermore,  in  the  earlier 
studies,  the  belt  forces  were  transmitted  through  the  shoulder  rigid  bodies  only 
to  Tl,  T2  and  T3,  further  increasing  for; <*>  these  vertebrae. 

In  the  initial  studies,  the  combined  belt  force  was  too  high,  because  more 
accurate  belt  force  time  history  data  were  not  yet  available.  Therefore,  rather 
than  trying  to  predict  injury  potentials  for  studies  1-4,  we  concerned  ourselves 
primarily  with  the  gross  effects  observed  in  these  studies. 

It  can  bo  seen  from  Table  6.2,  that  when  Tl  is  4 in.  forward  (PAM4),  the  neck 
axial  force  and  sagittal  plane  bending  moment  increased  about  40%  respectively  over 
the  2 in  forward  initial  configuration  (PAH2) , The  axial  forces  in  the  T3-T4  and 
L3-L4  discs  were  affected  vary  little.  The  sagittal  plane  bending  moments  at  these 
levels  increased  58X  and  33%,  respectively.  These  increases  in  sagittal  plane 
bending  moments  can  be  attributed  to  increased  inertial  affects  of  the  bead -helmet 
mass.  Comparison  of  the  150. msec  configurations  of  PAH2  and  PAH4  shww»  that  the 
upper  thoracic  region  of  PAH4  exhibits  considerably  more  sagittal  plane  curvature 


than  PAM2 . 
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The  results  of  Study  3,  which  considered  the  effects  of  varying  the  reel  height, 
are  presented  in  Figa  6.13,  and  6.14  and  Table  6.2.  Lowering  the  reel  height  from 
0.9  in  to  2.8  in  results  in  a considerable  increase  in  the  internal  axial  forces 
and  moments.  The  vertical  component  of  the  resultant  force  at  the  shoulder  rigid 
body  centroids  increased,  causing  the  crew  member  to  be  pushed  down  through  the 

. j . - .O  /'■  > . I f 

shoulder  rigid  bodies  as  well  as  back" into  the  seat. 

When  a portion  of  the  seatback  was  cut  off  so  that  the  head  and  vertebrhl 
bodies  T1-T5  do  not  contact  it*;  the  forces ^sustained  by  the  spine  increased  O 
dramatically.  If  we  compare  the  results  listed  in  Table  6.2,  we  find  that  the  peak  1 

i' 

axial  forces  In  the  neck  and  the  L3-rL4;disc  are  Increased  by  63%  and  86%  respec- 
tively, and  the  peak  sagittal  bending  moments  are  increased  by  406%  and,  20%  respec-, 
tively  when  the  seatback  is  ',eut-o£fM.  The  peak  mid-thoracic  axial  forces  and 
sagittal  bending  moments  are  increased  also.  The  100  msec  pilot  configuration. 

Fig.  6.15  corresponds  approximately  to  the  time  of  the  peak  neck  sagittal  bending 
moment.  In  this  configuration,  the  upper  thoracic  region  has  been  pulled  back  over 
the  top  of  the  seatback  while  the  head-helrntt  mass  is  still  slumped  forward. 

The  results  of  Study  5,  depicted  in  Fig.  6.17,  show  far  less  curvature,  of  the 
spine  because  of  the  reduced  belt  force  levels.  This  is  due  to  the  improved 
modeling  of  the  shoulder  rigid  bodies,  as  well  as  the  considerably  less  severe  belt* 
force  time  history.  If  we  check  the  internal  forces  against. the  injury  criteria 
presented  in  AMRL-TR-76-10,  it  is  found  that  internal  force  levels  in  the  spine 
(T1-L5)  are  well  below  injury  levels. 

The  inertial  effects  of  the  head-helaet  mass  are  particularly  evident, in  the 
150  and  200  msec  pilot  configuration  of  Fig.  6.19  from  Study  6.  Except  for  large 
forward  rotations  at  Tl,  the  spine  appears  to -behave  quite  well.  The  100  snd  150 
msec  configurations  appear  to  indicate  large  downward  rotation  of  the  rib  cage 
relative  to  the  spine*  due  to  the  increased  curvature  of  the  thoracic  region  of  the 


spine.  This  difficulty  arises-  because  there  is  no  obstruction  to  downward  rotation 

of  the  rib  cage  in  PAH  6,  the  inclusion  of  theabdominal  viscera  in  this  model 

would  limit  this  downward  motion*  ,t  ,,  v 

It  may  be  of  some  interest  to  note  that  the  15Q  and  200  msec  configurations  of 

studies  5 and  6 respectively,  are  very  similar  in  spinal  alignment  except  for  the 

larger  forward  rotations  at  Tl  in  Study  6.  In  Study  5,  the  belt  interacted  only 

with  the  shoulder  while  in  Study  6,  the  belt  also  interacted  with  rib  pairs  R1-R7.  t 
It  is  observed  that  in  both  cases,  the  spinal  configurations  when  the  seat  is 

contacted  are  similar.  From  Table  6.2,  it  can  be  seen  that  including  belt  inter- 
action with  the  ribs  has  the  effect  of  reducing  the  compressive  forces  in  the 

N 

lumbar  region,  and  although  the  sagittal  plane  bending  moment  Is  somewhat  increased 
in  the  lumbar  region,  it  appears  to  be  distributed  more  uniformly  along  the  spine. 

Figure  6.20  depicts  the  injury  potential  function,  IPF,  for  Study  6.  All  the 
vertebral  bodies  (particularly  in  the  lumbar  region)  are  well  below  the  injury  level, 
The  vertebral  body  Tl  has  not  been  shown  because  of  the  unrealistically  high 
neck  moment. 

In  summary; 

1)  The  forces  caused  by  pre-ejection  retraction  alone  are  not  sufficient  to 
lead  to  lower  vertebral  injuries;  however  it  does  induce  curvature  of  the  lower 
spine,  and  the  studies  of  Belytschko,  et  al.,  (1976)  have  shown  that  increased 
curvature  leads  to  increased  force  levels  in  ejection. 

2)  The  forces  and  initial  curvature  are  increased  as  the  reel  height  relative 
to  the  shoulder  is  decreased. 

3)  Isolated  spine  models  can  reasonably  well  reproduce  the  results  of  the 
more  complex  models. 

4)  The  forces  sod  curvatures  in  the  upper  spine  are  strongly  affected  by  tbs 

motion  of  the  neck  and  spine,  so  improved  models  of  the  neck  are  needed  to  obtain 
a complete  picture  of  the  response. 


• - CHAPTER  VII 

1 EJECTION  SIMUIATICNS 
VII. 1 Introduction  and  Obiectivea 

In  this  Chapter  we  consider  the  response  of  the  spine  models  is  ejection. 

These  studies  have  three  objectives: 

1)  to  evaluate  the  effects  of  various  components  of  tbe  model,  such  as  the 
viscera  and  seat  belt,  especially  with  respect  to  the  force  levels  in  th*  spine 
and  the  curvature} 

2)  to  compare  the  results  of  models  which  have  been  matched  with  experimental 
impedance  data  with  the  previous  models  in  Belytschko,  et  al.,  (1976) } 

3)  to  compare  the  results  obtainable  from  models  of  varying  degrees  of  com- 
plexity, such  as  the  simplified  spine  model,  SSM  and  the  complex  spine  model,  CSM, 
in  ejection  simulations. 

One  of  the  goals  of  the  second  objective  was  to  observe  the  effects  of  para- 
meters which  must  be  chosen  quite  arbitrarily  in  a simulation,  such  as  head 
positions.  It  should  be  noted  that  some  of  these  studies  were  made  before  our  final 
viscera  data  was  completed  through  the  impedance  studies, 

Except  where  noted,  the  acceleration  profile  used  in  all  of  the  studies  has  a 
rate  of  onset  of  714  g/sec  for  14  msec,  and  then  a constant  acceleration  of  10  g for 
66  msec,  for  a total  duration  of  80  msec.  After  80  msec,  the  acceleration  vanishes. 
The  acceleration  vector  is  perfectly  vertical  (4s)  for  ail  runs.  In  all  cases,  the 
problems  are  treated  3S  three  dimensional,  although  a two  dimensional  representation 

would  often  be  adequate.  The  explicit  integration  scheme  was  used  with  a time  step 
-4 

of  10  seconds,  thus  requiring  800  time  steps  for  a solution  duration  of  80  msec. 
VII  .2  Isolated  Ligamentous  Spine  Model 

VII. 2,1  Unrestrained  Ligamentous  Spine  Model 

The  first  two  simlatioue  involved  the  isolated  ligamentous  s*iae  models 
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(ILS  and  IlSV  from  Chapter  II)  and  considered  tha  effects  of  the  secondary  coluam, 
which  represents  the  ribs  in  the  thoracic  region  and  the  viscera  in  the  lumbar 
region,  cm  the  transmission  of  the  acceleration  induced  forces  along  the  spine. 

The  sect back  was  included,  but  the  restraint  system  was  not.  Table  7.1  gives  the 
peak  internal  forces  for  the  simulation  without  and  with  the  secondary  column; 

Fig.  7.1  depicts  the  initial  and  80  maec  pilot  configurations  for  the  simulation 
without  the  secondary  column  while  Fig.  7.2  depicts  the  80  msec  pilot  configuration 
with  the  secondary  column.  The  initial  pilot  configuration  was  the  same  as  in 
Fig.  7.1. 

From  Table  7.1  we  can  see  that  including  the  secondai  column  causes  a slight 
increase  in  the  axial  compressive  forces,  while  the  sagittal  plane  banding  moment 
is  reduced  considerably  in  the  lumbar  region  and  increased  slightly  in  tha  mid  to 
upper  thoracic  region.  The  secondary  column  is  a simplified  representation  of  the 
secondary  path  of  force  transmission  provided  by  the  viscera  and  ribs.  When  it  is 
absent,  the  spine  and  its  associated  ligaments  are  the  only  path  for  the  entire 
load.  Lucas  and  Bresler  (1961)  demonstrated  that  the  static  frontal  plane  buckling 
load  for  an  isolated  ligamentous  spine  constrained  against  sagittal  plane  motion 
is  2 x 10^  to  10  x 10^  dynes.  The  sagittal  plane  bending  stiffness  is  somewhat 
lower  chan  the  frontal  plane  bending  stiffness  due  to  the  sagittal  plane  curvatures. 
This  would  indicate  a lower  static  sagittal  plane  buckling  load.  Since  the  axial 
compressive  forces  transmitted  along  the  spine  during  an  ejection  simulation  are 
several  orders  of  magnitude  greater  than  the  frontal  plane  static  buckling  load  as 
determine  by  Lucas  and  Bresler  (1961)  one  would  aspect  a response  such  as 
depicted  in  Fig.  7.7  when  u*e  secondary  path  of  load  transmission  is  not  included. 

One  result  that  may  appear  contradictory  at  first  is  that  the  axial  load  is 
reduced  when  the  secondary  colusa  is  not  included.  However  in  the  absence  of  a 
secondary  column,  the  vertical  acceleration  of  the  head  and  upper  torso  are  reduced 
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Table  7.1  Comparison  of  peak  Internal  forces  for 
ejection  sit&ulationa  with  and  without 
secondary  coluan  (no  restraint  system^ 


Level 


JUial  force 

Q 

dyne 8 5 s 10 


Moment 

dyne -cm  x 10^ 


without 


without 


l-Kw^  ' 


-2.46 

-2.34 

-2.07 

-2.03 

-1.91 

-1.79 

-1.70 


-1.38 


-1.23 


-2.76 

-2.87 

-2.71 

-2.40 


-2.53 

-2.28 

-1.84 

-1.95 

-1.87 

-1.47 

-1.25 

-1.19 

-1.04 

-1.06 


-23.80 

20.60 

47.60 

57.87 

36.45 

-14.12 

-23.43 

-35.80 

-34.99 

-25.70 

-18.91 


-3.63 


-0.80 

-0.86 


-9.88 


'Refers  to  disc  below  designated  vertebral  level 


-6.48 

18.76 

28.18 

25.96 

22.56 

15.78 

-17.45 

-22.00 

-25.02 

-27.90 

-25.96 

-22.78 

-12.49 


10,21 

11.03 


considerably  In  the  fir.t  80  mane  taw.,  the  nxc.iv.ly  ««*u  spine  acts  as  a 
■'.hock"  absorber.  Thus,  whil.  the  aecondstycotancorrles  a signified  part  of 

the  toad  If  the  upper  torso  Is  constrained,  its  net  effo.fl.  to  Increase  the  anlet 
load  for  an  unrestrained  simulation. 

VII, 2. 2 Effects  of  Restraint  System  •: 

To  observe  the  effects  of  a harness  and  seatbelt  on  the  response  of  the  IIS 
model  (with  the  secondary  column),  an  ejection  was  simulated  with  the  harness  wodefcd 
as  tensile  springs  from  the  seatbsck  to  Tl,  T2  and  T3  and  the  seatbelt  as  a tensile 
spring  from  the  seatback  to  the  pelvis.  Ibis  simulation  was  rva  for  3200  time  steps 
for  a duration  of  320  msec.  Table  7.2  lists  the  peak  Internal  foreea  end  Figs.  7.3 
and  7.4  depict  the  0,  40,  60  . 80  and  120,  160,  240,  320  msec  pilot  configurations. 

As  expected,  the  restraint  system  haa  very  little  effect  on  the  axial  forces 
since  It  serves  primarily  to  redoes  lateral  motion.  The  sagittal  plana  banding 
moment  Is  almost  unchanged  In  the  lumbar  rerlo„,  but  Is  reduced  considerably  in 
the  thoracic  region  (e.g.  reductions  of  13%,  38%  and  SIX  at  T9,  T5  and  Tl, 

respectively  during  the  first  80  msec).  The  scatters,  as  modeled,  appears  to  have  no 
observable  effect  ou  the  response. 

It  is  noteworthy  that  some  of  the  peak  forces  occur  after  the  scat  accelerstlon 
goes  to  zero.  For  convenience,  wo  will  refer  to  the  first  80  msec  and  the  remaining 
240  msec  of  the  simulation  os  ^ and  AT,,  respectively.  boring  w2  the  axial  forces 
become  positive,  indicating  that  s tensile  w.  ,e  is  induced  by  the  seat  tccelorstion 
dropping  off  to  sere . The  peak  sagittal  place  bending  moment  magnitudes  at  levels 

T6  t0  11  ‘“'I  U,e  neck  occ'"  *»!»*  «r  Pot  T6  to  T3,  these  value,  are  an 
average  of  approximately  100X  gtestcr  than  too  respective  peak,  during  4T..  The  peak 
neck  moment  during  4I;  la  approximately  4.6  times  that  during  iTj.  If  „8  take  the 
nownt  in  Tl  to  be  given  approximately  by  the  overage  of  the  moments  in  the  T1-T2 
disc  and  the  neck,  this  would  yield  value,  of  4.92  x 107  and  16.98  , 107  dyuo-cm 


Level' 


Axial 
dynes  x 


Force 
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&>  7,2.  Feak  internal  forc^s  for  320  nsec  » jaot&ap^  ?«*• 
simulation  of  ILS  model  with  secondary  coluan 
; . and,,  restraint  system. : •*»'  •:  ••  «'•  Md?*#** 
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-2.29 
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-2.25 

18.27 

ni 

-2.52 

••>.  ,>  : .:•'!•  -17,26 

no 

-2.25 

,v.,,  ... ,^20.93 

T9 

-1.80 

-20.36 

T8 

-1.76 

-18.23 

T7 

-1.75 

-15.55 

T6 

-1.37 

5.70 

T5 

-1.29 

7.70 

T4 

-1.15 

6.60 

T3 

-0.83 

-3.72 

T2 

-1.17 

-5.43 

T1 

-1.15 

-3.64 

Neck 

-1.05 

6.20 

(l)Re£era 

to  disc  below  designs red  vertebral 

-9.82 

••  • . in  b 


12.61 

12.71 

12.81 

-8.55 

-6.94 

-5.49 
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Configurations  at  120,  160,  240  and  320  msec  for  ejection  of  ILS  with  restraint  system 
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Assessing  these  results  with  the  injury  criterion  ‘indicated  that  although  the 
combined  thoment  and  axial  force  Stresses  £or  tho?e  levels  discussed  above  arc'  be*ov' • 
the  failure  values  fdr  the  entire  duration,  at  Tl  the  possibility  of  injury  is  ' 
quite  high.  It  must  he  kept  in  mind  however,  that  in  this  mouel  the  neck  is  modeled 
as  a single  beam  element  between  the  head  and  Tl  and  considers  neither  the  effects  u- 
of  the  soft  tissues  nor  the  muscles,  which  would  most  likely  serve  to  reduce  the 
magnitude  of  the  neck  moment  transmitted  to  Tl.  The  modeling  of  the  harness  by 
springs’ 'connected  directly  to  Tl,  T2  and  T3  also  fails  to  account  for  the  much  less-  * 
severe  distribution  of  restraining  forces  through  the  shoulders  and  rib  cage  and 
associated  soft  tissues.  A larger  amount  of  damping  in  the  uodel  would  also  lower 
the  force  levels  during  This  simulation  demonstrates  that  the  introduction 

of  a restraining  system  decreases  the  flexural  response  of  the  Isolated  Ligamentous 
Spine  during  ejection,  and  permits  the  pilot  to  remain  upright.  The  larger  duration 
time  allowed  for  observance  of  longitudinal  and  flexural  frequencies  of  13  and  10  Hz 
respectively. 

VII. 2. 3 Increased  Lumbar  Stiffness 

Belytsehko  et  al.,  (1976)  noted  thi.t  if  the  preloaded  state  resulting  from  body 
weight  is  taken  to  be  the  reference  point,  a linear  approximation  to  a lumbar  disc 
load  deflection  curve  (above  the  preloaded  state)  results  in  a larger  stiffness 
than  from  the  unloaded  state  by  approximately  33%.  The  previous  simulations  did  not 
reflect  this  value,  so  a run  was  made  with  the  ILS  model  in  which  the  axial  and 
bending  stiffnesses  of  the  lumbar  intervertebral  discs  were  multiplied  by  1.333. 

The  simulation  was  the  same  as  that  previously  described  except  for  the  increased 
lumbar  disc  stiffness;  the  duration  of  the  simulation  is  60  msec. 

It  was  determined  that  the  change  in  the  axial  force  peak  magnitudes  is  essen- 
tially negligible  while  the  peak  sagittal  plane  bending  moment  magnitude  is  primarily 
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increased  in  the  lumbar  region  by  s*average<>£ 

the  spine  configurations  were  unchanged*  Thus  increasing  the  lumbar,  stiffness  by 
33%,  had,  negligible  effect,  sp  the  .spinal;,  .response,,  _ Paying,  thl^  ^raij.letiot^,,.t^a, 
maximum  combined  force  in  the  spring  elements  representing  the  ahoulder  hamess 

g 

was  approximately  1.65  x 10  dynes  (370  lbs).  . ., . ,• 

VII. 3 Ligamentous  Spine  with  Riba  and  Hydrodynamic  Visceral  Representation  . 

VII, 3,1  Normal  Configuration 

A 10  Gz  ejection  simulation  with  acceleration  profile  and  restraint  system  as 
described  in  Section  VII. 2 was  run  with  the  ligamentous  spine  with  riba  and  hydro-  ...... 

dynamic  visceral  representation  (CSM  from  Chapter  II),  Table  7.3  gives 

the  peak  internal  forces,  tUu  peak  visceral  pressures,  and  the  peak  harness  forces lP . . 

respectively.  Fig.  7.5  depicts  the  initial,  40,  60  and  80  msec  pilot  configurations 

respectively. 

This  simulation  was  exactly  the  same  as  the  one  in  Section  VII. 2 (for  the  first 
80  msec)  for  which  the  results  are  given  in  Table  7,2  and  Fig.  7.3.  Comparing 
Tables  7,2  and  7.3,  a small  increase  in  the  axial  force  is  noted.  The  peak  sagittal 
plane  bending  momenta  show  considerable  reductions.  In  the  lumbar  and  thoracic 
regions  and  the  neck,  they  are  reduced  by  averages  of  65%,  52%  and  21%  respectively. 
As  can  be  seen  from  the  oQ  and  80  msec  pilot  configurations.  Figs,  7.3  and  7,5, 
the  curvature  of  the  lumbar  region  is  less  for  the  hydrodynamic  visceral  represen- 
tation than  for  the  secondary  column  representation. 

5 2 

Table  7.3  indicates  a variation  in  peak  pressure  values  from  4,9b  x 10  byne/cm 

5 2 

in  the  lower  portion  (PELVIS -15)  of  the  viscera  model  to  2.51  x 10“'  dyne/cu  in  the 

upper  portion  (L3  to  the  bottom  of  rib  pair  10),  Morris*  et  al.,  (19bl),  using 

interabdominal  pressure  transducers,  determined  an  abdominal  pressure  of  2 at  IG3 
2 

dyue/ca  during  weight  lifting  experiments.  In  a dynamic  environment,  wecaa  expect 
response  magnitudes  of  up  to  twice  the  corresponding  static  pressure.  Hence, 
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Table  7.3  Peak  internal  forces  for  ejection 
simulation  of  the  CSM 

Level 

(1)  Axial  force 

4 : > ■"  O 

: . dvucc'x  10” 

Hcinent  9 

dyne-cmx  U/ 

J-  -J*  ' ■ U V ' ‘ * * 

L3 

,r  v-3.23  a \\ 

v.  V.  W-1.82  .. . 

. \ • -4  ’ 

14 

K\  "~\-  • 'y^  -3.00 

\ v \ ; ' T.85;\ 

A:-. 

L2 

' ~ ' v / 

«*..  f V •;  -K* 

^^.so---' 

Ll 

- '.wW 

-2.67>f^ 

'1  "J  JfLaV*.  ./•*'* 

X12 

-2.50 

/ 6.69 

Til 

-2,49  fs. 

-8.41 

XiO 

■"  V \ ■ 

'V  '.;v;  -9.93 

X9 

v.  ' >2.28  V 

>vr.\V*^‘i65 

18 

>,-2,39/  .,v  ’ 

-11.67 

17 

X6 

■ /<  .*  ;>^i.,95-;:  • 

';7>s  -7.-?9  }' 
.>>'y>4.12:"  " 

X5 

‘ 2.31 

X4 

-1.50 

3.05 

T3 

-1.03 

-2.95 

X2 

-1.39 

-2,46 

Tl 

-1.21 

• 1.21 

Neck 

2.12 

4.89 

Level 

• ' > : 

2 5 

Pres eore  (dvne /cm  x 10 

Lower (Pelvis  - 15) 

%>’’  -4.96 

Middle  (L5-L3) 

-4.06 

Upper 

(L3-bottoco  of 

-2.51 

rib  pair  10} 

Level 

P-orce  (dynes  x 10?) 

Pelvis  (seatbelt) 

ks7  ....V  . . 

T3  ■'] 
T2  j 

V (harness)  : 

3.10  ^ 

TW 

f 

1.51  / 

/I  V . 

v 'Refers  to  disc  below  designated  vertebral  level 
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Configurations  at  0,  40,  60  and  80  msec  for  ejection  of  CSM  (normal  configuration), 


comparison  ptorris^  et  al.  |^ul^  tn^a^^e^eak  visceral  pressure 

values  observed^  this  Wiinulatiom  ire  reasonable.  The  combined  peak  harness  force 
7 

is  9.52  x 10  dynes  (214  lbs)  in; this  simulation. 

VTI.3,2  Effect  of  Variation  in  Harness  Location 

In  the  previous  simulations,  the  harness  was  attached  to  the  seatback  at  the 
undeforaed  height  of  Tl.  To  see  Ihow  a change  in  reel  height  affects  the  response, 
an  ejection  was  simulated  with  the  initial  reel  height  at  the  undeformed  height  of 
T4,  thus  lowering  the  reel  height  by  6.3  cm  (2.5  in).  Table  7.4  gives  the  peak 
internal  forces,  the  peak  visceral  pressures,  and  the  peak  harness  forces 
respectively.  Figure  7.6  depiots  the  initial,  40,  60  and  80  msec  pilot  config- 
urations, respectively. 

Comparing  the  results  listed  in  Tables  7.3  and  7.4  indicates  that  lowering  the 
reel  height  from  Tl  to  T4  hac  the  following  effects:  the  axial  forces  in  the  lumbar 
and  thoracic  regions  are  increased  slightly  while  the  axial  force  in  the  neck  shows 
a small  reduction;  the  peak  sagittal  plane  bending  moments  undergo  neg^’gible  change 
in  the  lumbar  and  lower  thoracic  region  but  increase  an  average  of  701  from  T8-T4, 
decrease  by  4031  at  T3  and  T2,  and  increase  by  80%  and  100%  at  Tl  and  the  neck 
respectively. 

The  peak  visceral  pressures  and  the  total  peak  harness  force  show  very  little 
change,  although  towering  the  reel  did  change  the  distribution  of  the  harness  force. 
VII. 3. 3 Variation  in  Head/Helmet  Hass  Center  Location 

To  observe  the  effects  of  ehange  in  the  head-helmet  ©ass  center  on  the  internal 
force  response,  an  ejection  was  simulated  with  the  head-helmet  sags  center  saved 

forward  from  its  normal  location  by  4 in.  The  simulation  time  was  80  sec. 

Table  7.5  lists  the  peak  internal  forces  and  visceral  pressures.  Figure  7.7 

depicts  the  0,  40,  60  and  80  msec  pilot  configurations.  Figures  7.8  and  7.9  show 
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' Table  7,4  Peak  Internal  forces  for  ejection  simulation 


Level 

of  CSM  with  reel  located  at  the  initial  level 

of  T4,  " • v -•  - "'J 

Axicl  force  •'Moment  ' ■ 

S 9 

dyne  x 10  dvne-era  x 10 

15 

-3.31 

..  -1.67  . 

14 

-3.10 

7.81 

L3 

-3.07 

10.85 

L2 

-2.95 

10.02 

Ll 

-2.79 

8,58 

T12 

-2.69 

6.70 

Til 

-2.53 

-8.43 

TlO 

-2.54 

-10.08 

T9 

-2.42 

-8.65 

T8 

-2.70 

-13.01 

T7 

-2.34 

-9.43 

T6 

-2.10 

-6.36 

T5 

-1.76 

-6.48 

T4 

-1,74 

-3.96 

T3 

-1.34 

-1.51 

T2 

-1.44 

1.58 

Tl 

-1.19 

2.22 

Neck 

-1.04 

9.93 

Level 

2 S 

Pressure  (dyne/ca  x 10  ) 

Lower 

(Pelvis  - L5) 

-5.26 

Middle  (LS-U) 

-4.38 

Upper 

(LS-bottcaa  of 

-2.57 

rib  pair  10) 

Level 

7 

Restraieiag  torce  (dyaea  x 10  ) 

Pelvis  (seatbelt) 

9,69 

T3'4) 

12  ( 

(harness) 

0,57 

Tl-J 

3.30 

4,91 

^Refers  to  disc  below  designated  vertebral  level 
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Table  7.5  Peak  internal  forces  for  ejection  simulation 

of  CSM,  head/he lmet  mass  : an ter  moved  4"  forward 


Level 


<i> 


Axial  force 
dyneg  x 10 


Moment 


3 


yne-ca  x 10 


8 


L5 

' ■ -3.25  . . 

4,09 

L4 

-2.99 

7.98 

L3 

-3.00 

10.56 

12 

-2.88 

9.63 

Ll 

-2.62 

8.22 

Tl2 

-?.*6 

-7.42 

Til 

-2.47 

-8.51 

TlO 

-2.41 

-10.10 

T9 

-2.13 

-8.67 

T8 

-2.33 

-11.73 

T7 

-1.88 

-8.18 

T6 

-1.67 

-8.89 

T5 

-1.33 

-6.80 

T4 

-1.36 

-6.41 

T3 

-1.14 

-5.69 

T2 

-1.06 

-9.11 

Tl 

-0.88 

-10.00 

Neck 

-0.80 

-39.  n 

Level 

2 5 

Visceral  Pressure  (dyne7c»  x 10  ) 

Lower  (Pelvis  - L5) 

-5.02 

Middle 

(L5-L3) 

-4.56 

Upper 

(L3-bott«a  of 
rib  pair  10) 

-2.67 

^Refers  to  disc  below  designated  vertebral  level 
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Figure  7.8.  Axial  force  response  in  tfee  S-LS  «nd  tUtU  discs  and  the  neck  for 
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selected  axial  force  and  sagittal  plane  bending  mxaeat  titan  histories.  The  peak 
axial  compressive  forces  were  very  similar  to  those  observed  in  the  previous 
simulation  with  this  model.  There  was  some  decrease  at  Tl  because  the  head 
moved  forward,  so  Its  vertical -acceleration  was  reduced.  The  peak 
sagittal  plane  bending  moments  from  the  lumbar  to  the  mid-thoracic  region  era  also 
not  significantly  altered  from  previous  results.  In  the  upper  thoracic  region 
(T4-T1)  and  the  neck,  the  peak  moments  are  increased  considerably.  If  we  take  the 

axial  compressive  force  and  moment  in  Tl  to  be \the  averages  of  the  quantities  in  :~ 

7 8 • 

the  T1-T2  disc  and  the  neck,  8,4  x 10  dynes  and  2.5  x 10  dyne-cm  respectively, 

then  the  injury  criteria  of  Belytschko,  et  al#,  (1976)  indicates  that  the  stress 

level  in  Tl  is  well  above  the  failure  criterion.  The  80  msec  pilot  configuration. 

Fig.  7.7,  demonstrates  that  the  deformation  of  the  neck  is  quite  severe. 

The  peak  visceral  pressures  ir.  Table  7.5  are  little  changed  from  previous 
simulations.  This  is  expected  since  the  axial  forces  and  moments  in  the  lumbar  and 
lower  thoracic  region  also  showed  little  change. 

From  the  ejection  simulation  presented  in  this  and  the  previous  section,  we  may 
conclude  that  (1)  the  increase  of  the  curvature  of  the  spine  is  reduced  considerably 
by  the  introduction  of  a restraint  system  (i.e.  harness  and  seatbelt),  (2)  including 
a detailed  representation  of  the  rib  cage  and  viscera  also  increases  the  stability 
of  the  spine:  (3)  the  internal  forces,  particularly  the  sagittal  plane  bending 
moments,  are  effected  by  variations  in  inertial  reel  height;  (A)  a variation  in 
initial  haad/heltaec  mass  location  can  lead  to  substantial  increases  in  sagittal  plane 
bending  momenta  and  forward  motion  of  the  head;  (5)  the  axial  forces  predicted  by 
the  model  seem  quiet  independent  of  these  parameters. 


m- 


Wf 


m 


HR- 


1 


m 


ft- 


BT- 


m 


m: 


wmm 


m- 


wm 


StS-J'V 
Safes  • 


; 

• i 


Vll.4  Comparison  of  Ejection  Responae  of  SSM.  ILSV  and  CSH 

The  SSM  was  developed  So  reduce  the  computational  effort  of  various  simulations 
by  replacing  various  groups  of  elements  from  the  more  complex  models  by  single 
elements.  V©  have  already  shown  the  similarity  of  the  SSM  and  ILSV  impedance  curves. 
In  this  Section,  we  will  compare  the  SSM  response  In  ejection  to  the  ILSV  and  CSM 
models. 


An  ejection  simulation  was  run  with  the  SSM  using  the  standard  +10G  acceleration 

z 


profile  described  in  Section  VII. 1.  A restraint  system  was  included  consisting  of 
a seatbelt,  modeled  as  a tension  only  spring  between  the  seatback  and  the  pelvis, 
and  a harness,  modeled  as  a tension  only  spring  between  the  seatback  and  Tl.  The 
reel  height  was  the  initial  height  of  Tl.  The  ejection  environment  of  the  SSM  was 
therefore  identical  to  the  CSM  ejection  simulation  discussed  in  Section  VII. 3.1. 

A similar  ejection  simulation  was  run  with  the  ILSV.  The  initial  SSM  and  ILSV  con- 
figurations- are  the  normal  configurations  shown  in  Figures  2.7  and  2.4  respectively. 

The  damping  parameters  of  the  SSM  and  ILSV  were  adjusted  so  that  the  lowest 
spine  axial  and  bending  modes  and  the  lowest  visceral  mode  are  damped  by  10%.  This 
is  consistent  with  the  CSM  for  which  the  lowest  spine  axial  mode  and  lowest  visceral 
mode  are  also  damped  by  10%.  These  damping  levels  are  lower  than  those  determined 
in  the  impedance  study.  Justification  for  lowering  the  SSM  and  ILSV  damping 
values  is  based  oa  the  work  of  Vykukal  (1965).  Vykukal  exposed  four  subjects 
in  ? semi- supine  position  to  a vertical  acceleration  of  £ 0,46  in  a frequency  range 
from  2*}  to  20  Ha  combined  with  linear  accelerations  of  1,  2*  and  4G.  He  reported 
that  the  impedances  of  the  subjects  Indicated  a decrease  in  damping  with  an  increase 
in  sustained  acceleration. 

Tables  7.6  and  7.7  compare  the  peak  internal  force  responses  and  visceral 
pressures  of  the  SSM,  ILSV  and  CSM  during  the  standard  ejection  environment.  The 
peak  axial  force  In  an  SSM  element  is  compared  with  the  average  of  the  peak  axial 
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forces  in  the  corresponding  discs  of  the  ILSV  end  CSM,  The  peek  SSM  moments  ere 
taken  from  the  upper  end  of  the  element.  For  example,  the  peak  moment  for  element 
L3-T10  is  taken  from  the  T10  end.  This  is  compared  tc  the  peak  moment  at  the 
superior  end  of  the  specified  disc  of  the  ILSV  and  CSM.  The  peak  visceral  pressures 
of  the  SSM  visceral  elements  is  computed  by  dividing  the  peak  internal  force  of  the 
element  by  the  average  cross-sectional  area  of  the  element.  Hence,  this  is  actually 
the  peak  visceral  element  stress.  The  corresponding  value  for  the  ILSV  is  the 
average  of  the  individual  peak  stresses  in  the  noted  ILSV  visceral  elements. 

The  peak  axial  forces,  moments  and  visceral  pressures  show  good  overall  agreement 
both  in  magnitude  and  distribution.  Some  discrepancy  results  from  the  lumping  of  the 
masses  in  the  SSM.  Although  not  tabulated,  the  peak  harness  forces  show  reasonable 
agreement  for  the  three  models. 

The  peak  vertical  accelerations  at  Tl  are  17. lG,  20. 2G  and  18. 8G  for  the  SSM, 

ILSV  and  CSM  respectively  (G  «=  acceleration  due  to  gravity).  The  peak  vertical 
accelerations  at  the  head  are  20. OG,  24. 0G  and  19. 7G  for  the  SSM,  ILSV  and  CSM 
respectively.  The  agreement  is  quite  good. 

Although  not  shown,  the  40,  60  and  80  msec  pilot  configurations  for  the  three 
models  are  vary  similar. 

The  results  presented  in  this  section  demonstrate  that  the  SSM,  ILSV  and  CSM 
overall  ejection  simulation  responses  are  in  good  agreement.  Bence,  the  SSM  say 
replace  either  the  ILSV  or  the  CSH  for  ejection  studies  in  which  the  general  aspects 
of  the  response  are  of  interest  rather  than  the  detailed  res  pease  at  each  vertebral 
level.  The  SSH  would  be  vary  useful  for  design  studies  where  short  running  time  is 
a premium,  such  as  in  ai  interactive  mode. 
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This  Appendix  briefly  describes  the  programming  procedures  for  obtaining  the 
Impedance  curves  of  the  spine  models. 

AI.l.  Description  of  Procedure 

The  impedance  curves  of  the  spine  models  are  of  interest  in  the  frequency  range 
of  approximately  the  0-30  Hs.  This  is  well  below  the  numerical  stability  limit  for 
explicit  time  integration  of  the  spine  models.  For  this  reason  and  the  fact  that 
impedance  is  obtained  from  a linearized  analysis  (thus  requiring  only  one  assembly 
of  the  stiffness  matrix),  the  implicit  code  is  used.  The  procedure  for  obtaining  an 
impedance  curve  is  outlined  as  follows: 

(1)  in  subroutine  FREEFD,  a unit  step  velocity  is  prescribed  at  either  the 
pelvic  or  seat  primary  node  global  z translation; 

(2)  in  subroutine  SCLY1,  this  degree  of  freedom,  called  NS  EAT , must  be 
identified.  A double  precision  array,  FSEAT(1025)  is  included  to  provide  storage 
for  the  reaction  force  time  history.  At  the  end  of  each  solution  step,  the  statement 

FSEAT  (NSTEP)  - -FINT  (NSEAT) 

stores  the  negative  of  the  internal  axial  force  acting  on  the  pelvis  or  buttocks  in 
FSEAT.  The  acceleration  of  the  driven  mass  is  zero  except  at  time  zero  (the 
velocity  is  a step  function),  so  FSEAT  represents  the  reaction  forte; 

(3)  after  the  last  solution  step  in  S0LV1,  the  FFT  (Fast  Fourier  Transform)  sub- 
routine is  called  by  the  statement 

CALL  FFT<F3EAT , DSLT ,KS£STEP) 

where  BELT  is  the  ablution  time  step,  and  MXSTEP  is  the  number  of  solution  steps; 

(A)  subroutine  FFT  computes  the  DFT  (Discrete  Fourier  Transform)  of  FSEAT,  The 
modulus  of  impedance,  which  la  defined  as  2(a)  * iv(ttj) * determined  by 


dividing  the  magnitude  of  the  BFT  of  FSEAT  by  the  magnitude  of  the  DFT  of  the  unit 
step  velocity  at  each  frequency  step.  If  the  unit  Step  velocity  ia  applied  directly 
to  a mass,  (e.g.  at  the  pelvis)  the  scalar  PMASS  must  be  set  equal  to  the  trans~ 
lational  mass.  The  impedance  of  the  mass  is  then  also  computed.  If  the  unit  step 
velocity  is  not  applied  directly  to  a mass  (e.g.  at  the  seat),  PMASS  must  be  set 
equal  to  zero.  The  FFT  impedance  output  is  in  the  form  of  F(K) , Z(K),  ZTOT  and 
ZMASS  at  frequency  step  K, where 

F(K)  = circular  frequency  at  frequency  step  K, 

Z(K)  = modulus  of  impedance  (not  Including  the  impedance  of  the  driven  mass) 
at  frequency  step  K, 

ZTOT  *•  total  modulus  of  impedance  at  frequency  step  K • Z(K)  + ZMASS, 

ZMASS  = modulus  of  impedance  of  the  driven  mass  at  frequency  step  K. 

AI,2.  Data  Notes 

A method  for  computing  the  approximate  duration  time,  AT,  required  for  the 
transient  response  to  approach  zero  is  described  in  Chapter  IV.  Once  AT  is  deter- 
mined, the  accuracy  of  the  FFT  depends  on  the  number  of  solution  steps,  MXSTEP.  It 
was  found  that  MXSTEP  * 256  is  a satisfactory  number  of  solution  steps  for  obtaining 
impedances  or'  the  spine  models  in  the  0-30  Hz  region.  Hence,  the  solution  time  step, 

BELT,  is  given  by  AT/256.  The  upper  bound  to  MXSTEP  is  1025  since  all  arrays  are 

N 

dimensioned  at  this  value.  If  MXSTEP  # 2 , }J  £ 10,  a cubic  interpolation  routine, 

subroutine  INTRP,  realigns  the  discrete  points  in  the  time  history  so  that  it  is 
N 

defined  by  2 points,  N s 10.  If  passible,  however,  it  is  recommended  that 
N 

MXSTEP  * 2 since  this  will  eliminate  the  interpolations. 

For  the  impedance  runs  as  well  as  for  other  runs  in  which  damping,  greater  than 
approximately  101  of  critical  is  employed,  it  ia  recommended  that  the  stiffnos* 
proportional  damping  is  used  rather  than  the  fraction  of  critics*  damping  because 
this  damping  behaves  better  for  larger  damping  values. 
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